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Problems of flow, heat release, and mixing in gas turbine combustion chambers are dis- 
sussed and analyzed. The analysis is based on an idealized or equivalent chamber in which 
the end result is the same as in the conventional chamber, but in which the functions are 
separate and distinct A large primary zone is advantageous, but the combustion process 
still is burdened by an inverse relationship between primary air and fuel. Conditions favor- 
the primary zone over the whole range of operation are realized by 
means of a Mfiechanical control Agreement between analysis and experiment was demon- 
strated by tests of a chamber with a mechanical control. Comparison of the field of mixing 
of primary flame gas and secondary air with other studies of mixing revealed a general simi- 
larity, but a slower rate of mixing that is ascribed to the effects of a large difference of temper- 


able to combustion in 


ature or density in these experiments 


1. Introduction .= Length of potential core, ft 
Absolute static pressure, lb/ft 

. . Absolute total pressure, lb/ft? 

h jet propulsion for aircraft has become a } V2: dynamic pressure, Ib/ft® 


reality, many problems connected with it Radius, ft 
ily partially solved. One of these problems Radius at which T (or V) is halfway between T 
elopment of combustion chambers for turbo- _(or V) at axis and that of secondary stream, ft 
. ‘ Gas constant, ft-lb/lb ° I 
s that will give satisfactory performance T= Absolute temperature, ° R 
service conditions. Satisfactory perform- ‘= Velocity, ft/sec 
es that a continuous supply of exhaust gas Weight flow per unit time, lb/sec 
esired temperature is furnished the turbine, t= Axial distance from station 3, ft 
i. An fw ofl ion oft enadll y = Lateral distance, ft 
response to demands for changes of condl- z= Upper limit of probability integral 
is positive and occurs In a minimum of time. Emissivity 
lesirable characteristics are high combustion Combustion efficiency 
ey, small loss of pressure, uniformity of tem- Density, slugs/ft 
tor he e | a oe id icht Scale factor 
ire in the exhaust, minimum size and weigh enteient shone ainees Moitt 
with required performance and life, and 
of manufacture and maintenance. Some com- 
ses have been necessary in designing combus- 
chambers, since some of these features have been Fuel 
1 . Indicated 
nd to be obtained only at the expense of others. Dlnsaey 
Problems of flow and combustion in conventional Secondary 
bustion chambers are so complicated and inter- 
ted that development of chambers has progressed 3. Combustion Chamber 
by empirical methods. Although develop- 
t by these methods has produced immediate re- 
is, in the long run a fundamental understanding | developed and modified according to the require- 
the problem is advantageous and will promote | ments of the system in which they are intended to 
It is the purpose of | be used. Although each chamber may be different 
from any other, it appears that most designs differ 
in detail only, and not in principle. An examina- 
sin an attempt to arrive at a better understanding | tion of the conventional combustion chamber reveals 
their influence on the performance of the chamber | that there are four main functions that it must 
perform, (1) separation of the air into primary and 
2N 1 secondary streams, (2) maintenance of a stable pilot 
- Nomenciature zone, (3) combustion of fuel with the primary air, and 
(4) mixing of the burned gases with secondary air. 
Area, ft? \ "eke, : , . 
E@ectivences of mixins Design features may be illustrated by the schematic 
Coefficient in mixing length equation drawing in figure 1. A protected, turbulent zone of 
Specific heat at constant pressure, Btu/Ib°F. low velocity ts established at A inside of the dome of 
Diameter, ft : the perforated liner, and when fuel is injected into 
Standard acceleration of gravity, ft/sec?. . 
Mixing length, ft this zone a mixture is formed which is capable of 
" | sustaining a stable flame. More air enters the liner 


ribed in this report was sponsored by the Bureau of Aeronautics 
e Navy es between A and B and combustion proceeds as it 
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Combustion chambers for gas turbines have been 
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s paper to separate and examine some of the func- 
ns of conventional gas turbine combustion cham- 
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mixes with fuel. Velocity of the mixture in this zone 
Is probably higher than the velocity of the flame in 
the mixture, but combustion is piloted by the flame 
from A. The air that takes part in the combustion 
in this zone is known as primary air. Secondary air 
enters and mixes with the products of combustion 
between B and C, thereby reducing the gas tem- 
perature to a safe operating level. Completeness 
of combustion is promoted if the air is admitted 
gradually into the liner and in a manner so that it 
does not quench the flame and cause blowout. 

There is no assurance that all of the air that enters 
the liner between A and B takes part in the combus- 
tion process, and therefore it may be said that the 
separation into primary and secondary air is not 
distinct 

Gas turbine combustors must be capable of oper- 
ating over a wide range of exhaust to inlet temper- 
ature ratio, and generally this flexibility is attained 
by injecting atomized liquid fuel into the pilot zone 
at A This method of fuel injection prov ides 
heterogeneous mixtures that are capable of sup- 
porting stable combustion, even though their quality 
may differ substantially from stoichiometric. It 
does not follow that combustion must always be 
efficient with this system. Actually, it may be good 
over only a narrow range of conditions, and may 
become poor when others are imposed. 

Loss of pressure due to flow friction is high, and 
problems of flow and combustion are complicated 
and interrelated. Flow is turbulent, and in some 
regions within the primary zone, it may actually be 
in an upstream direction. Combustion and mixing 
within the liner nowhere occur as distinct, individual 
processes. Although these conditions appear to be 
beneficial from the standpoint of performance and 
flexibility, they make a detailed analysis of the 
conventional chamber almost an impossibility. 

If a chamber is designed to carry out the functions 
enumerated above as separate and distinct opera- 
tions, then an analysis of its performance becomes 
possible. A chamber designed for this purpose may 
be called an idealized or equivalent chamber, 
equivalent in the sense that the end result is the 
same as in the conventional chamber. The dif- 
ferences between this chamber, represented by 
figure 2, and one of conventional design are apparent. 
Stations A, B, and C in each figure represent locations 
where the same respective operations are performed. 
In this equivalent chamber the air is separated into 
two distinct streams, and the drop in pressure due to 
flow friction that each experiences is made to occur 
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in a manner that is subject*to analysis. A thot: 
process is assumed for the primary air, and 
secondary air is assumed to go through a throttliy 
process followed by a reversible expansion. Fy 
injected and combustion is stablized after th prima 
throttling process, and is assumed to be comple; 
within the primary zone. The size of the pilot yo, 
has been diminished to that which would exist ; 
flame holder, and combustion may be said to oery 
as inramjet chambers. Since combustion is made 
occur at a fixed temperature ratio in a stream | 
has a uniform, measureable velocity, the process 
subject to treatment by both analytical and expe 
mental methods. Mixing of the burned gas wy 
secondary air occurs only after comb Istior 
complete, and the changes resulting from this pro 
also may be predicted and measured 


4. Analysis of the Equivalent Combustio; 
hamber 


The method used in the analysis depends on t! 
fact that each function of the combustion chambe 
is considered to be performed as a distinct, individu 
operation. Changes of temperature, pressure, et 
are calculated for each, and from these the perform 
ance of the chamber as a whole is derived 

Before these calculations may be made, it is neces 
sary to specify some of the conditions at which th 
chamber is required to operate. Although an 
practical chamber must be able to operate over 
wide range of conditions, there is usually some con 
dition where it performs best, and it is desirable ths 
this should coincide with the condition that is used 
most often in service. For the purpose of this ana 
ysis, an exhaust to inlet temperature ratio of 3.28 
a suitable specification, because it is close to that 
used in practice and is also a convenient working 
condition for experimental purposes. 

As combustion is usually best when the mixture ¢ 
air and fuel is at or near stoichiometric, this condition 
may be used as the basis of a second specification 
A mixture of this quality would give an exhaust 
inlet temperature ratio of about 7.25 across | 
primary zone. 

Some loss of pressure in conventional chamber 
occurs when the air flows through the perforation 
of the liner. In order to approximate this condita 
in the chamber represented in figure 2, it is necessary 
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loss of pressure in the flow of air in a 
its effects can be determined analyti- 
xperimentally. Once the magnitude of 
pressure due to friction is specified, the 
mn of the chamber may be calculated to 
f the above specifications, provided the 
if the flow in the chamber are also speci- 
ever, the solution of the problem would 
to the set of conditions of flow that were 
Because an infinite number of conditions 
and as all of them cannot possibly be 
is necessary to make a choice. The 
simplified if it is that the 
flow in the chamber approaches zero, and 
ption is used throughout these analyses. 
yn chambers of gas turbines are operated 
vy low velocities; hence such an assumption 
cause serious errors 
mplifying assumptions may be made with- 
serror. These are: rise of temperature in 
wry zone may be accomplished without 
fuel to the flowing stream; specific heat, c,, 
is is independent of temperature; loss of 
due to friction does not occur except where 
elocity of flow is uniform in each zone across 
the duct normal to the flow; and secondary 
ves no heat before it mixes with the gas at 
perature from the primary zone. 
ie problems of flow, heat release, and mixing 
be handled by using Bernoulli’s equation and 
servation of momentum equation. Since it 
msidered that the velocity of flow approaches 
Bernoulli’s equation for incompressible, fric- 
ess flow may be written 


considered 


re p, and p are constant at all points in the flow 
conservation of momentum equation for heat 
ise OF Mixing 1s 


Pr pV? constant, 


when heat release occurs where V approaches 


eq 2 may be written 


subseripts 1 and 2 in eq 2a indicate the state of 
gas before and after heat release, respectively. 


When these equations are applied to cases where V, 
is appreciable, it is equivalent to considering that the 
fluid is incompressible, or that density is dependent 


only on temperature. 


The condition for continuity 


of flow may be written 


Eq: 
suff 
the 


pAV 


constant. 


lations 1, 2, and 3 and the gas Jaw, p=gpRT, are 
icient for the analysis of the problems of flow in 
equivalent combustion chamber. 


4.1. Effect of Size of the Primary Zone 


Two design conditions of temperature rise in the 
chamber have been specified, but nothing has been 
said about the size of the primary zone. It would 
be possible to specify arbitrarily a size based on 
results of practical experience, but there is no assur- 
ance that the decision on this basis would be the 
most logical A determination of the effect of size 
of the primary zone by analytical means probably 
would prov ide a sound basis on which to make a 
choice of size. This determination is attempted in 
the analysis that follows 

A list of the design the 
numerical subscripts refer to the stations in the 
chamber indicated in figure 2, is presented below. 
The subscripts p and s refer to primary and secondary 


conditions, in which 


air, respectively. 
(1) 2,/7,;=3.2: (2) 73.,/7,;=7.25; (3 


, APs/Q2.9 
2=loss of pressure due to flow friction; (4 


2 AP »/Qe.p 
10=loss of pressure due to flow friction; (5) ps,.=pas.»- 

At low velocity, loss of pressure due to friction 
may be expressed as a constant number of velocity 
heads, g, where g=4pV*. The losses indicated in 
the third and fourth specifications were chosen in 
an attempt to approximate those in practical cham- 
bers. The fifth specification is that pressure is 
uniform across the plane of station 3. Tollmien [1] ° 
calculated that the pressure at the plane of the 
discharge of a free, round jet differs from the am- 
bient by an amount equal to about one-tenth per- 
cent of the dynamic pressure of the jet. A pressure 
difference of this small magnitude would have no 
significance in the present analysis. A variation of 
serious magnitude probably would not be encoun- 
tered except where the velocity of flow in the primary 
stream at station 3 equaled the acoustic velocity 
in the stream. Conditions such as this are not 
encountered in present-day engines, and consequently 
are not considered here. 

The above specifications provide sufficient infor- 
mation to permit a calculation of the effect of liner 
size, d., on the problems of flow, heat release, and 
mixing in the equivalent chamber. 

If W, and W, are the rates of flow of the primary 
and secondary streams, respectively, then when 
mixing of these streams occurs between stations 3 
and 4, a heat balance may be written 


W,(73,.— T,)e,>= W,(TM.— Te. (4) 


J 


Noting that Cy 1S assumed to be independent of 
temperature, eq 4 may be written 


and when the first and second specifications are 
substituted, the right member of eq 5 is equal to 


——————- 


2 Figures in brackets indicate the literature references at the of this paper 
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Fiacure 3 


0.35. In other words, 35 percent of the air must 
enter the primary zone to fulfill the conditions of 
specifications (1) and (2 The air that enters the 
primary zone can be regarded as the central core 
of the stream at station 1. Since it is considered 
that the properties of the flow are uniform at this 
station, the diameter of this central core of primary 
air can be expressed by the equation 


W, 
4/DY= We Ww, 


Bernoulli's equation, eq 1, may be used to deter- 
mine the changes of pressure between stations 1 and 
2 as d, is changed. Thus 


V3 
Pi— P2= 42-9 1 (3-1) 


and the ratio of the velocities for both primary and 

secondary air may be determined by the continuity 
equation as follows 

D\* 

p) 

d? TP 2 

hal (9) 

d? TP ‘ 


Since (d,/D)*’=W,/(W,+W,), the pressure changes 
may be expressed by the equations 
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These were evaluated and plotted versus (d,/)) 
figure 3,a. When (d,/D)* is greater than 0.35, th 
primary air diffuses into the primary zon 


| velocity decreases, and its pressure increases. Whe 


d,/D)*=0.6, the velocity of the primary air as estal 
lished by eq 8 is about 58 percent of V;. 

Loss of pressure in the primary zone is the sum 
that due to friction and combustion. This may | 
written 

(p2— Ps)p = Ap, +242 »( uf -— ] ); 
l, 
where the second term on the right side is determine 
by eq 2a, which is the loss due to combustion on th 
basis of the momentum equation. 4q,,/q, is given bY 
eq 8, and the pressure loss thus may be expressed 1 
terms of g, by the following equation: 


(Po—Ps)p_[Ap,, o(Ta» di/p\' 
-P : [ S242 r |G °) 


P 





valuated and the results are plotted in | of course, on the relative pressure losses in the pri- 

mary and secondary zones. For example, if no 

m of the change of pressure in the | friction loss were assumed, the maximum area would 

air is based on the assumption that | be 66 percent, whereas if both the assumed losses 

In order to fulfill this condition, it is were doubled, the maximum area would decrease to 
o change d; with d,, and the required | 58 percent of the total 

y be evaluated by using eq 1 and 3. At station 3 the primary air at high temperature 

pressure in the secondary air may be | and the secondary air begin to mix, and the resulting 

changes of pressure may be calculated by the equa- 

14 tion for the conservation of momentum, For this 
case it may be written 


iv be solved for p Eq 7 and 12 may be d 
0 yield prt+2e0(1—4) +0424 


Pe Pp q ' APs, 


1), 15) | This equation may be arranged to the more con- 
venient form 


is considered that ps, ’z,., the right P Ps_o Ws 9 @ d _ 9 13.3 d,\’ 
. , ini = (1-35) “ (7) 


15 may be substituted for ps,, in q qi qi 
this method the following relation is 


er 


and the change of pressure is determined by evalu- 

ating each of the separate terms and taking the 

Pi— Ps) + (4 Ap,) —(pi-— 2) »- 16) | sum. These terms may be evaluated by using the 

continuity equation and the gas law, and always 

may be related to gq, by the continuity | considering that p is dependent only on temperature. 
and this relation leads to Thus 

» Ws 9 Tr ( T, Dp ) » rs D ), (20 

); 17) en TAT, &) ~? TAG, 20) 


9 13.3 d? - D?— d?\ d? 
D)*, the following is * ( l= 7: ) 2( d—d: ) ( l—7 


j p Ia0(% ‘ 9 i? (“ T ' (4 y ” T ; d‘ . 
. ~ a \d a Ge TF d ~ dy GG 
Dp: j d? 
( 4 ) ig The values of these terms were determined from the 
q previous relations, and the sum expressed as (p3—p, 

q, is plotted versus (d,/D)? in figure 3, d. The over- 

all change of pressure from station 1 to 4, which was 

evaluated by a summation of the successive changes, 

is plotted versus (d,/D))? in figure 3, e. 


D?— d? 
d? 


PUNE BS 


(3) 


Geta ieriwes se we 


ese equations furnish sufficient information to 
the change of pressure and of d; as d, is 


Thus eq 14 and 16 evaluate (p..,—ps; 
nd eq 16 and 18 evaluate (d;/D)))*.. Changes of pres- Although the over-all change of pressure is informa- 


between stations 2 and 3 are plotted in figure tive, it does not allow ready comparison with the 
and (d;/D)? is plotted in figure 3, ¢. Both are | change of pressure of other combustion chambers 
d against (d./D)*, which defines the size of the In order to make a comparison on this basis, it would 
mary zone Figure 3, b, shows that the pressure be necessary for each chamber to have inlet and 
ss decreases rapidly as the size of the primary zone exhaust ducts of the same relative size However, 
creased. However. the velocity in the secondary | if the loss of pressure is expressed in the fundamen- 
becomes high, and a point ‘is reached (when | tally more useful form of total pressure, a comparison 
0.6) where the loss in pressure across the | '5 possible without this information about size. Loss 
is so high that there is no need for more drop | of total pressure may be written 
essure by expansion in the nozzle at station 3. 
spointd;=D. If the primary zone were made Ps_ 444 4 
ull larger, the loss of pressure in the secondary air q qi 
vould be so high that an enlargement or diffuser 
Thine ould be necessary to recover and thus maintain the | The first two terms on the right 
On eRuality of pressure across station 3. Thus d; would | previously (see fig. 3, e, and eq 20), and the loss of 
ven "¥ have to become larger than D, and this is shown in | total pressure was determined and is plotted versus 
ssed }, ¢, by the curve of ()/d;)*, where this ratio | (d,/D)* in figure 3, f. Inspection of this curve 
weomes zero when (d,/D))?=0.65. For all practical | reveals the advantage of a liner of large size. When 
irposes, it may be said that the area of the primary | the liner of this chamber occupies 60 percent of the 
this chamber cannot exceed 60 percent of the | total area, the loss of pressure is about half that when 
a. This area of the primary zone depends, | the liner occupies 35 percent of the area 
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side were eV aluated 





Little is known about the factors that control the 
rate of mixing of gas streams under conditions that 
exist in combustion chambers. In practical chambers 
the rate is accelerated by the turbulent flow of the 
air through the perforations of the liner. Although 
this is acc omplished at the expense of loss of pressure, 
it is compensated for by the consequent reduction in 
length of the chamber. One of the factors that may — le 
influence the rate of mixing is the ratio of the velocity : i 
of the gas streams at station 3. The ratio of q3 ,/q3,; Ficune 4. Diagram of mizing feld 
may be determined from eq 21 and 22, and noting 
that 7 T,=7.25, the velocity ratio may be derived, 











of primary to secondary stream velocity 
, decreased and approached unity. Beyond / 

( . )- d ) of decay of velocity along the axis was not ; 
d,\?\D greatly by V;,,/V:,, and therefore for the p; 


and equals 


/ d present purposes L may be regarded as ar 

—( dD ) measure of the effectiveness of the mixing process 
bringing the two streams to a uniform stat, 
other words, at a given value of axial distanes 


This was evaluated and the ratio is plotted versus 
- . x >L), decay along the axis would be less complet 


(d,/D) in figure 3, g. It increases with an increase of cca te telial 

“ ° en L 1s large. 
the size of the liner, notwithstanding the fact that ei Tr yee 
V,,, and therefore V » decreases The degree of > — om rouncer present values of L 

2, Bs & Veolt and in this analysis the variatior 
expansion required in the secondary stream also ; k 

; V,,/V;,. with d./D has been determined. <A com)j. 
decreases, and V,, decreases at a greater rate than : ; 
V nation of the two sets of information w 
determination of Z as a function of d., an 


; : presented in figure 3, h. In their analysis the constant 
incompressible mixing of a round jet in a general : , 
(’ was regarded as independent of J  y J and this 


stream, and the results were used here in an attempt ; 
‘ , , : » | same assumption is made here. At small values of 
to determine the effect of liner size on the mixing of 
d is also small, and mixing would probabl 
the primary and secondary streams. Although the od - : 
: : . good As d, is increased, V; ,/ V3, increases slow 
conditions existing in a combustion chamber differ 4 
- | but not sufficiently to counterbalance increase of 
in many respects from those assumed in the analysis, : 
.~. ’ | and Z increases. However, a point is reached wl 
it is possible that the two cases would exhibit similar > , — 
wt ar vee V; ,/ 1 increases at a rate sufficient to counte: 
Ip ; ; : balance increase of d,., and ZL begins to decreas 
In the analysis [2] the momentum transfer theory - : . = ' 
: <= | When (d,/D)?=0.6, L is well below the maxim 
was used to determine the way in which velocity is , 
om E ; m4. f and it may be said that mixing probably would 
transferred between two round streams Turbulent Fe 5d : . 
relatively effective at this condition. In connect 


shear stress, 7, defined by the equation 
with these conclusions, it must be remembered t! 
dV\? there are differences between the conditions assume 
pl ( ) 2: in the theory and those existing in a chamber 
di large difference of temperature, confinement of | 
mixing region by the walls of the chamber, and 
was equated to the change of momentum flux of the | relatively large primary stream would proba 
mixing region with axial distance. Mixing length, /, | influence the mixing in a manner that would req 
was assumed equal to the width of the mixing region | experimental determination 
multiplied by a constant of proportionality, C. The This analysis has indicated the effect of the siz 
constant ( was considered a characteristic of the | the primary zone on the problems of flow, heat 
turbulent motion, and appeared in the end result as lease, and mixing in the equivalent chamber 
a constant to be evaluated experimentally. It was | primary zone of large size appears to offer the advai 
also assumed that the variation of velocity across | tages of low velocity in the primary zone, and 
the mixing region was similar to that determined by | duced loss of pressure across the chamber. Low 
the cosine equation velocity probably would simplify the problem of 
maintenance of a stable flame, and a reduction in t! 
loss of pressure would increase the ‘‘mechan 
efficiency” of the chamber and thereby increas 
over-all efficiency of the engine in which it would | 
where r,, is the radius at which velocity is half way | used. With regard to the problem of mixing, th 
between axial and secondary stream velocity. Re- | analysis indicates no serious disadvantage attached 
gions in the flow and significance of terms in eq 26 | to a large primary zone. Maximum area of th 
are illustrated in figure 4 primary zone in the equivalent chamber is found t 
Extent of the potential core of undisturbed pri- | be 60 percent of the total chamber area. Measure- 
mary stream, L/d,, was found to increase as the ratio | ments of practical chambers indicate that the ares 
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ry zone is generally about 60 percent of 
d this point of similarity between analy- 
tice 18 gratifying For these reasons the 
primary zone is fixed at the maximum 
analysis that follows 


alysis of Operation of the Equivalent 
Combustion Chamber 


means of control over the over-all tem- 
tio in conventional chambers is by the 
f the rate of fuel injection into the pri- 
This may be approximated for pur- 
alysis by assuming that the temperature 
T;, in the primary zone is allowed to 
hile holding the configuration of the cham- 
int, as determined Before pro- 
» the calculations, it is advantageous to list 
ng complete specifications for this combus- 


above 


ber 

0.6D?; (2) d,=D; (3 ADs/q2 2—loss of 

lue to flow friction; (4) Ap,/q2,,>=10=loss 

ire due to flow friction; (5) p3.=ps.»; 

varies with T 33 T;. 

lem is to evaluate the distribution of air 

the primary and secondary zones on the 
hat ps will always equal p Exhaust to 
mperature ratio, 7,/7,, can be determined 

distribution is known. Loss of pressure in 
mary zone is given by eq 12, and by using 

may be written in terms of p,, as follows: 


P2.7 Pp q q 
ng these two equations gives 


(2? T° 


Ps) q rT <4 T -1)—4. 28) 


f pressure in the secondary air is given by eq 
{and by evaluating p,, in terms of p, and com- 
ng, there results 


) 


Pi~— P. <4: q q 


ng 28 and 29 on the assumption that p 


12.5 (APs +2 I, ‘ n 30) 
q2; 7 


is specified equal to D, q..=4qs,,, and eq 30 
da relation between g,», d2,, and 73,,/7), as 


q23 (Ap, dos) 


lation, when combined with the equation of 
tv of flow, 


+(p VA), , (32) 


will yield a relation that will predict how the air 
divides between the primary and secondary zones as 
T;,,/7T, is changed. Thus the continuity equation 
may be written in the form 


and since V, ,/V»,,=(q2.s/q2.p)'”, the fraction of 


going into the primary zone is 


14 Ap, 72,97 . 
Ap 


Equation 5 expresses the relation between 7,/7; and 
W,/(W,+W,), and this may be combined with eq 
34 to yield 


This expresses the relation between the overall 
temperature ratio and the temperature ratio in the 
primary zone, and is plotted in figure 5,a. The 
relation between 7/7; and W, (W, W,) is plotted 
in figufre 5, b 

Inspection of these curves reveals some interesting 
information concerning the operation of the chamber 
Although 7; ,/7; increases with 7;,/7,, the fraction of 
air going into the primary zone decreases. This may 
be interpreted to mean that as more fuel is injected 
into the primary zone, the flow of air into it will 
decrease. Thus the length of the flame will increase 
with 7\/7;, and at values higher than the design 
condition it may extend downstream beyond the 
primary zone. Mixing with secondary air may not 
be complete by the time the gas enters the turbine, 
and the turbine may be exposed to gas at excessively 
high temperatures. A further possibility is that the 
conditions may not be favorable for combustion 
beyond the primary zone. The flame may be 
quenched by the large quantity of excess air, and 
combustion efficiency will then be low. In addition, 
the pilot zone may become so rich in fuel that it will 
not support combustion, and rich blowout will occur 

Conversely, at low values of 7,/7,, the mixture in 
the primary zone will be lean in fuel, and combustion 
may be poor. As 7,/7; is decreased, a point will be 
reached where the mixture in the primary zone will 
be too lean to support combustion, and lean blowout 
will occur. Moreover, the velocity in the primary 
zone increases and adds to the difficulty. These 
effects are observed in testing of practical chambers, 
although not exclusively because of the reasons 
enumerated 
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FIGURE Effects 


In contrast to the air distribution determined 
above, a distribution is shown based on maintenance 
of near stoichiometric conditions in the primary 
zone. Equation 5 was used with a constant value of 
T;,,/T; equal to 7.25, and the dotted curve tabeled 
“ideal” in figure 5, b,is the result. This ideal distri- 
bution is radically different and points out the fact 
that the conventional chamber has proper distribu- 
tion only at the design condition, that is, at the value 
of T, 7; where the two curves cross 

Both the type of mixture of air and fuel and the 
velocity of the air in the combustion chamber in- 
fluence the lean and rich limits of combustion. 
Heterogeneous mixtures are used exclusively because 
the range of combustible limits is relatively wide 
compared to that of homogeneous mixtures. It is 
yrobably true that low velocity tends to widen the 
imits in all chambers of this type 

When the distribution of the air is known, the 
changes of pressure and velocity in the chamber may 
be determined. Because the fraction of air going 
into the primary zone changes with 7/7), it is evi- 
dent that the streamlines of this air will diverge or 
converge on going from station 1 to 2, depending on 
the quantity that is flowing. If this fraction is 
smaller than that of the fraction of the area included 
by the primary zone, the streamlines will diverge, and 
the air will diffuse into the primary zone. Change of 
pressure between stations | and 2 may be determined 
as before by eq 5 to 11, remembering in this case that 
the fraction of primary air is the variable. A plot of 
these is presented in figure 5, c, for both primary and 
secondary air streams. At all values of 7,/7, the 
primary air in this chamber diffuses into the primary 
zone, whereas the secondary air expands. 


tT 


in the chambe 


with a large primary zone 


The changes of pressure between stations 2 a 
are evaluated by combining eq 28 and 10 for 
primary zone, and 28 and 11 for the secondary 
Thus by this method 


and 


These changes were evaluated, and the results 
plotted in figure 5, d. For this chamber the loss 
static pressure from station 2 to 3 does not chang 
greatly with 7/7). 

Mixing of the primary and 
begins at station 3, and the change of pressur 
results may be calculated as before by eq 19 throug 
22. This is plotted in figure 5, e, and the 
change of pressure from station 1 to 4 is plotte 
figure 5, f. Loss of total pressure may 
evaluated as before (see eq 23), and is plotted 
figure 5, g. In this equivalent chamber the 
loss is about 64. 

Observed loss of total pressure in a practical co 
bustion chamber is also shown for comparis 
These data were observed when V, was equal to 7 
ft/sec, and displacement of the curves from each ot 
is an indication of the relative magnitude of the 
of pressure due to friction. Loss of pressure 
combustion chamber is undesirable, because it low 
the over-all efficiency of the engine in which 
used, and from this standpoint it appears that t 
equivalent chamber would be superior. Howeve! 
size and performance of the chambers in other ! 
spects must also be considered before an effect 


secondary str 


ALSO 





over-a 


an be made. Combustion efficiency, 
are very important, but all the factors 
e these are not thoroughly understood 
as pressure, temperature, character of 
of air and fuel, etc. may not have the 
n all combustion chambers, and these 
idered in this analysis. Their influence 
hamber is determined almost exclusively 
nt 
of the velocities of the primary to sec- 
at station 3, V; ,/V was evaluated as 
eq 24) and is plotted versus 7,/7, in 
It increases with 7,/7;, and probably the 
ess will be most effective in bringing the 
niform state when outlet temperature is 
tive mixing is needed most at this con- 
oh temperature 
bustion chamber described and analyzed 
to practical chambers in which over-all 
mperature Ze T, is controlled by rate of 
on into the primary zone. In the analysis 
«1 of control is simulated by variation of 
re ratio, 7;,/7;, in the primary zone. It 
winted out that this system of control has 
desirable effects in the combustion chamber 
conditions of operation deviate from the 
int, the quality of the mixture of air and 
primary zone also deviates from the qual- 
s most favorable for combustion. Hence 
cy of combustion probably would decrease, 
eme conditions the mixture would become 
rich or lean to support combustion, It is 
that this system has certain advantages, 
simplicity of control, because all that is re- 
a control of the rate of fuel injection into 
ry zone. However, the advantages of 
ng conditions for good combustion are also 
and some effort may well be made to 
the requirements for this objective 


Analysis of Operation of a Chamber with a 


Controlled Primary Zone 


ted before, combustion is usually best when 
ire of air and fuel is near stoichiometric, 
e 7 7 may be selected to be equal to 
before The problem then is to select a 
“i and determine the conditions that are nec- 
to hold 7; ,/7; constant while 7,/7; is varied 
the first method that comes to mind is the 
of the size of the mixing section, since this 
was used in the first analysis. In that anal- 
as varied, and d,; was changed simultane- 
hold 7 5 7, constant This method would 
n adjustable nozzle at station 3, and in serv- 
bably would be impractical. However, an 
le throttle, or valve, also may be inserted in 
ndary zone for this purpose. Although this 
f control introduces loss of pressure due to 
t has certain advantages from the stand- 
mechanical simplicity. A combustion cham- 
such a device is represented by figure 6, 
throttle takes the form of two perforated, 
The pressure drop in the secondary 














= ADJUSTABLE COMBUSTION 


THROTTLE 


SECTION X-x 
Diagram of the chamber with the controlled primary 


zone 


FiGgure 6 


stream js adjusted by rotating one annulus relative 
to the other, thus providing larger or smaller open- 
ings as desired. This is only one of many possible 
schemes that might accomplish the same object. It 
will, however, serve for the purposes of this analysis 

The specifications for this case are (1) 7), ,/T 

7.25: 2 de 0 35D: 3 d DD: } Ap, Q: 
Ps. s=Ps.p; (6) T/T, varies with Ap,/q. 

Stations in the chamber are shown in figure 6 
It will be noted that d, was specified small (0.35'D)), 
and this was done for several reasons. Experimental 
work was planned on this chamber, and it appeared 
that conditions in a small primary zone would be 
more favorable in the proposed experiments to a 
well-reculated combustion process In a large zone 
the velocity of the primary stream would be low, 
and probably the flame would seat at the fuel in 
jector under some conditions, and under others 
would require a flame holder for stabilization. In a 
small zone, velocity would be relatively high, and 
combustion could be well regulated and performed 
in a homogeneous mixture of air and fuel with the 
seat of the flame at any desired position A large 
primary zone would require a throttle in the pri- 
mary stream, with attendant increase of mechanical 


problems and unknown effects on the flow condi- 
Another reason for 


tions and combustion process 
a small zone was that a study of the mixing process 
of flame gas and secondary air could be made under 
conditions more nearly comparable with those as- 
sumed in theoretical studies ot spread of free jets 
and of mixing of coaxial streams. The particular 
size chosen for the primary zone is that for which 
the flow of air will be undisturbed when operating 
at the design condition of 7,/7 A At this de- 
sign condition primary air is 35 percent of the total 
see eq 5 With regard to the specification of no 
pressure loss due to friction in the primary zone, 
this is not a serious departure from reality because 
the loss can be made small when combustion is per- 
formed under the conditions mentioned. These 
specifications of the method of combustion and of 


287 


Cha Mee cess 


arse are 
ed | 


we 


Ween? mer ewes 5 









































7 Effects in the chambe / the controlled prima 


pressure loss in the primary zone, although practical | across station 3, ps, 
in the laboratory, probably would not be practical | determined as follows 
in actual engines where flexibility of operation is 
important Pp Ps)s Pp Pp 
The fraction of air going into the primary zone 
under the conditions specified may be calculated | All of the terms on the right side have been evalua 
again by eq 5, using 7; ,/7, constant at 7.25. Re- | and their sum is plotted in figure 7, ¢ 
sults of these calculations are shown in figure 7, a The magnitude of this pressure change is var 
Change of pressure in the primary and secondary | over the operating range, and the variation is a 
air between stations 1 and 2 may be calculated by | plished by changing the opening of the adjust: 
using eq 10 and 11, and these changes are shown in | throttle shown in figure 6. The pressure chang 
figure 7, b expressed in terms of qg,, but it is also of interes 
Loss of pressure in the primary zone is define the change of pressure with respect to ¢ 
this may be of aid in the design of the thro 
die , Equation 9 is the relation between gq, and q, ,, al 
may be applied for this purpose as follows 
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and in order to evaluate 38 it is necessary to express Ap, _(p ; \P2 Ms (d./D ] 
72,» mM terms of q By Bernoulli’s equation, eq 1, Jos , d/D 
may be written 
This relation was evaluated and (p s) «/92 
plotted versus 7/7, in figure 7, d. When | 
throttle is opened, the fraction of air going into U 
and since (p,—p.), is determined by eq 10, eq 3 primary zone decreases, thereby decreasing 7 
and 39 may be combined to vield If the throttle were removed entirely (Ap 
, Z¢ T, would be about 1.8 Lower values ol f 
tl 2 could only be achieved by decreasing T;.; T, of 
( iv iw ) insertion of a throttle in the primary zone. 
axa . Change of pressure due to mixing of the secondar 
(: *) and primary streams may be calculated by eq 
D to 22. The term d; in these equations for this cas 
- is equal to D, and when this substitution ts ma 
his was evaluated and is plotted in figure 7, c. (pPs—p)/q: can be evaluated. Figure 7, e, 1° 
As it is considered that the pressure is uniform | plot of this quantity. The over-all change of pre 
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which was evaluated by adding the 
hnanges, 1s plotted in figure 7, f 
tal pressure may be determined by eq 23 
expressed by eq 20, is substituted in 23 
¢ relation is derived 


Ap 
q 


valuated and is plotted versus 7,/7; in 
Inspection of this curve reveals that the 
stantially higher than in the previous 
higher values of 7,/7;. The previous 
ileulated on the basis that pressure was 
» friction in the primary zone, but in this 
was assumed. In any practical applica- 
iamber of this type, some means of stabi- 
flame would be necessary, and this would 
loss of pressure in the primary zone. The 
of the loss of pressure required in service 
vould be comparable to that in the con- 
hamber. If this loss had been considered, 
ill loss of total pressure would have been 
in indicated in figure 7, g, and higher than 
conventional chamber 
o of the velocities of the 
reams at station 3, V;,,/V3,,, was evaluated 
otted in figure 7, h. This ratio 
/. and if this increase is a sign of an increase 
ctiveness of mixing, uniformity of temper- 
the exhaust will be good in this chamber 
nditions of maximum temperature rise 
roblem of control in combustion chambers 
tvpe would be more difficult than that of a 
mal chamber. A mechanism would be re- 
» adjust the throttle in the secondary air in 
to change of 7, or to change of over-all air- 
) (Another method, analagous to that used 
ocating engines, would be manual control of 
‘ttle, with fuel rate automatically controlled 
rdance with throttle opening. The problems 
facture, installation, and maintenance would 
e complicated, but improved performance over 
range might compensate for the added 
al complications 


primary and sec- 


increases 


5. Experimental Program 


specifications for the preceding analysis were 
with the thought that experimental work 
« undertaken to check the results of the 
It is relatively simple to match experiment 
ilvsis insofar as the physical size of the cham- 
temperature rise therein are concerned, but 
tion in other respects probably would influence 
its. The mere fact that combustion occurs 
addition of mass (of fuel) to the flowing 
Specific heat is not independent of temper- 

oss of pressure due to friction cannot be 
and exchange of heat between primary and 

ry streams before mixing cannot be pre- 
Although uniform velocity can be provided 
cross the plane of the chamber entrance in the ex- 


periments, the flow will take up a nonuniform dis- 
tribution at other stations. Operation of the cham- 
ber at a finite velocity of flow would also introduce 
some deviation due to compressibility, but not of 
serious magnitude if the flow Mach number is not 
Influence of compressibility on the flow 
some 


above 0.1 
relations is discussed in detail in reference [3] 
deviation between theory and experiment is to be 
expected, but probably is not of serious magnitude 


5.1. Experimental Apparatus 


The experimental combustion chamber was con- 
structed of standard 6 in. pipe with a 4-ft coaxial 
liner, 3%.— in. inside diameter and a wall thickness 
of %sin. The pipe was cut and flanged as illustrated 
in figure 8 for accessibility to the liner. A schematic 
drawing of the system, in which the location of the 
stations is shown, is also presented in figure 8 
Dimension X was 10 in. in the experiments on flow 
and heat release, and 8 ft in the experiments on 
mixing. Uniform velocity across the plane of the 
entrance to the chamber was provided by a nozzle 
of 12-in. diameter at the inlet and 6 in. at the outlet, 
or chamber, end. Throttles were placed in the 
secondary or primary zones as required to control 
the air distribution between zones. The throttle 
in the secondary zone was made of two thin annulli 
containing equally spaced %-in. holes. Rotation 
of one annular plate with respect. to the other con- 
trolled the opening and thereby the distribution 
Throttles in the primary zone consisted of disks with 
various sized holes at their centers 

Flow of air to the chamber was metered by a 
calibrated orifice located in the system far upstream 
from the burner. Liquid propane flowed from a 
weighing tank through a rotameter to a vaporizer 
tube in the shape of a hairpin mounted in the system 
far downstream from the liner. Gaseous propane 
from the vaporizer flowed to the multihole fuel 
injectol located in the liner about 6 in. downstream 


Fieure 8&8 The experimental chamber 
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from station 2. The injector was made of three In the experiments in which mixing of p 
concentric but not coplanar rings of tubing mani- | secondary streams was studied, most of 
folded together and drilled with numerous small | instrumentation was removed to prevent j 
holes. The injector shown in figure 8 was used | with the flow distribution in the cham 
only in the preliminary experiments and was made | siderable trouble with asymmetrical flow 
of two concentric rings. Mixing of the fuel with | in the mixing zone was experienced, a 
primary air took place between the injector and | necessary to rebuild the system on th 
flame holder, and the flame was initiated by a spark | The wall thickness of the liner was increa 
at the flame holder. The spark was turned off | in., and the liner was mounted coaxially 
after smooth combustion was established and the | length of 6-in. pipe Dimension \, fi 
flame was stabilized at the flame holder. The flame | these experiments was 8 ft. Elimination 
holder was a ring of triangular cross section, 2%-in. | joints in the vicinity of the liner by this m: 
outside diameter and 1-in. inside diameter. It | substantial improvement in the mixing x 
was placed 18 in. upstream from the end of the | liner of heavier construction warped less 
liner (station 3 A water-cooled valve was placed | also probably helped to improve the con: 
at the outlet of the pipe system, and this was used | the mixing study. 
to control the chamber pressure Instrumentation in the mixing zone of th 
Initial attempts to operate the chamber were | consisted of water-cooled total-pressure t 
uniformly unsuccessful because of rough or inter- | Chromel-Alumel thermocouples of 22-gag: 
mittent combustion. Numerous changes in the | iridium—50 percent iridium, 50 percent ri 
primary zone and downstream portions of the system | (Ir-IrRh) thermocouples of about 20-gag: 
made no improvement whatsoever, and it was not | Three of the Chromel-Alumel thermocouples 
until an orifice plate was installed 5 ft upstream from | driven by electric motors to traverse the pip: 
the nozzle that smooth combustion was possible. | speed of 2 in./min. This was equal to th 
This orihes plate was repla ‘ ad by one perforated speed of the temperature recorder, and by this t 
with numerous smaller holes, and all experiments | a full-scale profile of temperature was obt; 
reported herein were made with this plate in the | One motorized couple was located 32 in 
system. It is probable that this improvement was | 48 in. from station 3 and at right angles to ea 
due to change of the effective length of the column | Other Chromel-Alumel thermocouples we 
of gas in the system operated, and were located in the mixing 
Instrumentation was added for the measurements | tween the motorized couples 
of temperature, pressure, and flow in the chamber In order to measure the higher temperatur: 
Total-pressure tubes were installed at stations 1 and | mixing zone, Ir-IrRh thermocouples [4] 
2p. A pitot-static tube was also used at station 2p. | cooled supports were used. One of thes 
Static-pressure wall taps were located at stations | couples with its support tube is shown in 
0,1, 2p, 2s,3s,3p,and4. Temperature was measured | In this instance the wires are bent at right: 
at station 0 and also at a point far downstream from | the tube in an effort to gain depth of imm: 
station 4 by silver-shielded Chromel-Alumel thermo- | the gas. Only one bent couple was used 
couples. The pipe system contained two right-angle | projected straight into the stream at least 
bends beyond station 4 to insure that the gas was | These Ir-IrRh couples were operated by han 
well mixed by the time it reached the last thermo- | were used in the first 32 in. of the mixing zon 
couples The reading of this thermox ouple was con- 
sidered the average temperature of the exhaust gas, 5.2. Experimental Procedure and Results 
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Ve Some tests were conducted with insulation 
around the pipe beyond station 3 to determine the 
magnitude of the heat loss from the system. Tem- 


The burner was ignited at an air velocity 
tion 1 of about 20 ft/sec, and after smooth con 
tion was established, the rates of air and f 
increased simultaneously until the velocity at st 
1 was about 60 ft/sec. Measured temperat 
station 4 was maintained nearly constant during | 
operation because this would indicate “arly 
stant fuel-air ratio in the primary zone (at ag 
throttle opening Injection of the propan 
through the multihole injector provided n 
homogeneous mixtures of air and fuel. Lin 
smooth combustion in homogeneous mixtures 
narrow, and the method of control was usefu 
keeping within the limits while changing condi! 
in the burner. After the burner was brought 
the desired operating condition, final adjustme 
were made of fuel and air rates and burner pressu! 
level. The fuel was adjusted until the burne 
sounded smoothest. This criterion did not ne 
Ficure 9. /Jridium—iridium-rhodium thermocouple. | sarily guarantee the same fuel-air ratio in the primar 


peratures in the system were high, and the pipe was 
not insulated in the majority of tests because of 
possible damage to the system by the prevailing 
high temperature 
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sts, but it was adopted because opera- 
se (under conditions of resonant com- 
ld affect adversely the accuracy of the 
ts of flow and Although 
is not eliminated completely, the data 
its effects were reduced to a tolerable 


pressures 


he primary zone was operated a. nearly 
<ture ratio, the desired range of condi- 
er-all ‘temperature attained 
ise of throttles 

f air flow 25 |Ib/see in the 
Temperature of the air varied 
to 610° R, and pressure generally was 
itmospheres absolute Some variation 
onditions was tried but with no effect on 


rise was 
was about 1 


tests 


inlet to the chamber 

conditions of flow was about 1.2 in. of 
mperature ratio in the primary zone 
s calculated to be in the range of 6.6 to 7.1, 
es of 0.06 part of gaseous propane and | 
by weight over the range of values of 7 


q at the 


pressure 


work on the chamber was divided 
one in which problems of flow and 
studied and the other in which 

primary and secondary streams was 
Results of the former will be described 


‘ ntal 
arts, 
se were 


Experiments in Flow and Heat Release 


tioned before, fuel rate was adjusted in 
eriments until combustion was smoothest, 


equently the fuel-air ratio in the primary 


not constant in all tests. Tests in which 


the metered fuel-air ratio was 0.06 + .005 
lected, and the results of these tests are reported 
herein Selection by this rule reduced the spread of 
the data somewhat and did not 
relations 

Values of measured average temperature (7) after 
mixing of the gases by passage around two bends in 
the system indicated either low combustion efficiency 
S85 to 90 percent or large heat Analvses 
of the mixed outlet gas in three tests by a gravimetric 
method [5] indicated a combustion efficiency of 98 
percent, while a heat balance based on (i and the 
assumption of an adiabatic system indicated effi- 
ciency of about 90 percent. When insulation was 
installed around the pipe beyond station 3, efficiency 
based on 7, ranged between 92 and 100 percent 
and it was concluded that loss of heat from the system 
when uninsulated large. Prevailing tempera- 
tures in the system were high, and it appeared advis- 
able to operate without insulation for 
safety. On the assumption that the low values of 
efficiency were caused by heat loss and not by poor 
combustion, data in this report are based on tempera- 
ture T; calculated on the L00-percent 
combustion of the mixture of and fuel 
Cool secondary air contacted the walls of the 6-in 
pipe down to and somewhat beyond station 3, and 
it is expected that loss of heat from the gas down to 
station 3 and somewhat beyond would have 
negligible The majority of the observations were 
made in this upper section of the burner, and this is 
considered a justification for ignoring the heat loss 


were se- 


change the basic 


losses 


was 


reasons of 


basis of 


total als 


been 


Results of the tests are presented graphically in 
figure 10, and results of the last analysis (fig. 7) are 
included for comparison. The portion of the ai 
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Results in the experimental chamber 
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that entered the primary zone, figure 10, a, was 
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determined from pressure q@,, given by a pitot-static | values of (p,;—p.), were always correct 


tube at station 2p in the primary zone, or by the | chemical analyses also indicated essential 


fit 


difference of pressure from a total head tube and | plete combustion (based on oxygen cor 
static pressure wall tap. The data are displaced | across station 3p, and it was concluded 
from the theoretical curve, but compare well with | flame holder was not too close to station 3 
curves based on known combustion characteristics Change of pressure between stations 2 and 


¥ 


of propane-air mixtures. These curves were cal- | secondary zone is plotted in figures 10,¢,and 10 
culated for fuel-air mixtures of 0.055 and 0.065 in | the data. although displaced from the theor 
the primary zone, which is the range in the tests | curves, are in good agreement with the theor 


, 
z 
4 

~ 


reported These curves tak into account the sures between these two stations were also rr 

variation of specific heat with temperature and the | in the primary zone, and these are compared t 

added mass of fuel predicted by the conservation of momentum equ 
Division of the air into primary and secondary | in figure 12. In this figure observed va 


aéry 


streams between stations 1 and 2 is accompanied P2— Ps) /G2, » Were plotted versus fuel-air ratio 

by pressure changes, and results of measurements of | primary zone, and the theoretical curve is base 
these are presented in figure 10, b. The data are | calculated values of flame temperatures of prop 
again displaced from the theoretical curve, but the | air mixtures. Velocity, and therefore dynami 
same general relationship as predicted by Bernoulli's | sure, in the primary zone is increased by the addit 
equamon is evident Some of the displacement is | of fuel to the air stream. and observed Ja.» Was 
for reasons enumerated above. A comparison of the | rected as indicated in the figure to account for t 
measured data for the primary stream with the curve | Observed pressures were not corrected for th 
calculated from eq 1 is presented in figure 11, where | mentum of the fuel, which in these experiments was 
the line represents the equation In general it | injected downstream with a pressure drop in the 
appears that the equation is applicable. Because | jector of from 20 to 60 lb/in2. With no burning, 
magnitudes of the pressures are small, and any | served pressure drop was about 2¢,,,, which is du 
effect of resonant combustion in the chamber on | friction in the primary zone. This value of 29 

the measurements would be large on a percentage | added to the theoretical to account for the frict 
basis, it is difficult to decide which of the readings | drop. In the neighborhood of fuel-air ratio of 0.0 
are correct. Several chemical analyses of the exhaust | observed pressure drop is in good agreement with t! 
gas at station 3p were made, and the fuel-air ratios | theoretical. but at lean condition the observed value 
calculated from these analyses were in good agree- | are low. The reason for the low observed values ® 
ment with those calculated from values of q, derived | not clear, but it could be either incomplete combu: 
from eq | and the observed value of (p,;—p2),/q:. This | tion in the primary zone or an effect of resonance. If 
is not regarded as conclusive evidence that measured | combustion were not complete in the primary 20! 
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completion between stations 3 and 4 be- 
ed values of 7; normal com- 
ency 

bast assumptions in the analysis of the 
combustion chambet that pressure 
across the plane ot station 3 Some de- 
} this observed it 
yressure of the secondary stream exceeded 


indi ated 


was 


condition was was 


primary stream by about 0.1 q3,, in most 
rv as much as 0.3 q However 
the static pressure wall tap may hav: 


In some 


d adversely by warping and scale forma- 
by the high temperature of the liner. In 
appears that the pressure differences are 
to detract from the value of the 
basis of equal pressure. <A 
primary jet is indi- 
in the secondary 


nough 
the 


secondary air into the 


strong 
he surplus of pressure 


f the two streams takes place beyond sta- 
d the 
and 4 are plotted in figure 10,e 
the theoretical indicates that the observed 


changes between 
Compari- 


observ ed pre ssure 


hanges are large probably because of wall 
ic turbulence caused by the velocity differ- 
veen the primary and secondary streams 
vent of the friction drop with no burning 
it to be about 25 percent of the dynamic 
of the stream, which is a normal loss for this 
}in. pipe. During burning, average condi- 
flow across the mixing section would vive a 
pressure equal to g,7\/7\, and hence the 
on drop on this basis would be equal Lo 
ls \q If the data at the condition 7,/7,=4 
is an example, the observed (p3;—p,)/q, was 
ind correction by the frie tion drop A Pp q 
l’;=1) would lead to (p3— p,)/q,=0, which is 
ipared to the theoretical value of 2 The 
»is probably due to high level of turbulence 


king streams 
} 


in the experimental cham- 
and compares well with 


ll loss of pressuré 
tted in figure 10 


ated values 


b. Experiments in Mixing 


ombustion chamber was operated over a 

values of 7/7, of from 1.5 to about 4 by 

the size of the throttle and its position 
onding values of the velocity ratio, V3 ,/V3 >, 
> to about 0.1. The mixture in the primary 
is adjusted until it was on the lean side of 
metric by observation of the flame tempera- 
ile decreasing the fuel rate from the condi- 
Observations of tempera- 
lean 


reac- 


f excessive richness 
vere made in the mixing zone after the 
tion was established to reduce effects of 
the flame gas with secondary air during 
Only one measuring instrument was used 
1e to avoid disturbance of the flow. 
racy of the temperature measurements of the 
es was not of a high order, because only a 
alibration of the Ir-IrRh thermocouple was 


because corrections to the indicated 
temperature were large. At high temperature, heat 
by radiation is large, and consequently the 
junction would be at a temperature lower than that 
Corrections were calculated by a method 


available, and 
lost 


of the gas 
that is based on the assumption that the heat lost 
from the thermocouple by radiation is equal to that 
forced The 


corrections were cle l¢ rmined Is 


cained by convection equation by 


which the 


indicated bv the thermo 
couple, 7... was assumed equal to the g tempera 
ture next to the the pipe h is Stelans 
radiation constant (1.7*107° Btu/ft?/°R* hn = 
is the coefficient of transier Kor 
bare thermocouples, h, Btu/ft?/°R ho 
which is an empirical relation to take into account 
the effect of mass rate of flow G expressed in pou ids 
\ quick determination 
IrRh thermocouples 
Chromel-Alumel 


l’; is the temperature 


wall of 


convective heat 


225yG6 a) 


per square foot per second 
of the emissivity, e, of the Ir 
indicated it to be 0.26, and e for the 
the rmocouples was considered to be 0.9 

(As an example of temperature arrived at by this 
means, temperature of the flame gas at station 3 
fx. mm 13 averaged 4,350 °R with a 
mean square deviation of 140 ik This average 
temperature is about 300° F greater than the cal- 
culated flame temperature of a and air 
mixture near condition. It is not 
known whether the calculated corrections were 
large or whether the calibration of the thermocouple 
Some of the devia- 


tests root 


propane 
stoichiometri 


too 


is responsible for the difference 
tion mentioned was probably caused by inability to 
maintain the same fuel-air ratio in all tests 

Some measurements of dynami pressure were also 
made in the mixing zone, and velocity was calculated 
from these observed temper: 
tures and static pressures 

In the description of the experimental results that 
follows, the data with the ratio of 
velocity of the secondary to primary stream V {[F 
The values of this ratio reported were estimated from 
the observed rate of fuel injection (on the assumption 
of a constant fuel-air ratio of 0.062 in the primary 
and the observed temperature and pressure in 
The difference between these estimates 
than 10 
were 


measurements and 


were correlated 


zone 
the chamber 
and measurements 
percent Because the velocity 
not considered accurate, the estimated 
accepted as adequate for the present purpose 

1. Arial Distribution of Te m perature Mixing of 
the primary flame gas and secondary air takes place 
in a system bounded by the walls of the 6-in. pipe, 
and the gas will approach an average or mixed tem- 
perature 7 as it travels downstream from station 3 
Temperature at the axis of the burner, 7 was 
expressed in the dimensionless form (7 Ze 
7; ,—T, and correlated with distance z from station 
3 by an empirical equation similar to that derived by 


some was not larger 


measurements 
values were 
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Pai [6] for two-dimensional jet mixing. The equation 
involved the probability integral, as follows 


fi . 
’ d 5 
1 
The first member of eq 45 represents experimental 
observations , and values of 
were s« lec ted to make the two members of the equa- 


at various values of z 


tion equal 

In these experiments il Was not possible to make 
accurate measurements of 7,, because both velocity 
and temperature gradients existed in the stream, and 
no bends were used in the system to mix the gas 
Previous experiments in which mixed temperature, 
7,, was measured, indicated relatively large heat loss 
from the system These measurements were made 
downstream of the mixing zone, and if these were 
used as a guide to establish 7, in the present experi- 
ments, 7, by this method would probably be lower 
than the true value at station 4. Calculation of 7; 
on the assumption of 100-percent combustion effi- 
ciency, and no heat would give high values 
Both methods were used to establish values of 7;,, 
and although some difference was noted, the correla- 
tion obtained was not basically changed The wall 
of the pipe was relatively cool in the first half of the 
mixing section, and it seemed logical to conclude that 
T,, based on 100-percent combustion and no heat 
was nearer the true value, and this value of T; 
was used in the correlation to be described 

Values of 2 in eq 45 were determined from the 
observed axial temperature at various values of z, 
and these values of z were plotted against the corre- 
sponding values of z. It was found that plots of log 

against x/r;.,, Where rs, is the radius of the liner, 
could be represented by straight lines. Plots of this 
kind are shown in figure 13 for three different condi- 


loss 


loss 


20- 


The three lines ha 
and therefor 


tions of V; ,/V3.,<.1. 
mon point z=1.5, 
equation can be written 


rivs, y, 


where & is a constant, and ¢ is the slope o 
When 1; /,3.,>1, the same relationsh 
served, but only two out of five tests rave 
passed through the common point. 

A Maclaurin expansion of the probabili 
would indicate that (7, T)/(T3..—T, 
first approximation. This approximatior 
when 2<0.7, and in this range temperatu 
approximated by 


and from this it appears that ¢ is a measur 
effectiveness of the mixing in bringing the t 
ture at the axis to the mixed temperature 

Slope c was correlated with V3,,/V3., in figur 
and the data suggest a parabolic relationship, w 
a minimum at a velocity ratio of unity 

These results have been presented thus far wit} 
reference to or comparison with the theory of mix 
of jets. The mixing of streams, different in b 
temperature and velocity, has received very lim 
attention, and consequently the present results 
be compared with the theory for incompressibli fl 
presented in reference [2] and utilized in a previ 
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Correlation of axial temperature with distance 


25 


by the probability integral. 
T/Ti=23 





ld ft 


his report. A comparison will be pre- 
with the recently presented experimental 
Forstall and Shapiro [7] obtained in the 
axial streams at a common temperature 
d that differences shown by the compari- 
be caused by the large difference of Tem- 
density between the streams in the pres- 
ents. Hereinafter, comparison with the 
uire and Trouncer [2] is referred to as a 
with the theory, and the results of For- 
shapiro are referred to as the “‘other experi 


} is a comparison of the characteristics of 
as observed along the axis of the burner 
eory and also with the other experiments 
jet in a general stream would mix and 
approach the properties of the stream, the 
less fractions of temperature and velocity 
he temperature (and velocity ) of the sec- 
ream in order to make the comparison 
il, it may be said that in the median 
mperature expressed this way varies in- 
ith axial distance, and this is in agreement 
other experiments but not with the theory 
condition of | V5.» =0.375, limited infor- 
ibout axial velocity in the mixing zone 
that it is transferred at a slower rate than 
Three other tests at the condition 
also indicated this. Comparison of 
on velocity with the other experimental 
ndicates a somewhat slower rate of transfer 
esent experiments 
V3../Vs3.p was 0.125, the data show that 
of temperature in the present experiments 


ure 
t) 95 


mates that of velocity in the other experi- 


However, in the expectation that velocity 
transferred in this case also at a slower rate 


temperature, it appears that the rate of transfer 
ess in the present experiments than in the other 


rime 


nts 
the condition V;,/V; P 


One set of data points in this figure 
0.18, and this was 
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Ficure 14. Effectiveness of mixing. 


included because it illustrates more clearly the 
inverse relationship of temperature and distance 
In this system temperature approaches 7, as mixing 
becomes complete, not 7), and the completely mixed 
condition would be attained when (7 r 
(T: z T T,-—T Tr ; r The latter frac- 
tion is 0.47 when V3,,/\ 0.125, and the z data 
points show the beginning of a slower approach of 
temperature (at the lower values) to the mixed con- 
dition. At V;./V3.,=0.18, the fraction was 0.33, 
and the slower approach to the mixed condition 
occurred downstream from the place where the 
observations were made in this case 

It will noted that the straight 
extended to a value of unity for the temperature 
fraction, and at both conditions of V3,,/\/3,,<1 the 
axial distance is about llr,,. This same value of 
ll was noted when | V;, was 0.55, which 
leads to the conclusion that the burner can_ be 
assigned a characteristic length. The general inverse 
relationship of temperature and axial distance, and 
the characteristic length, would be of general utility 
in making an approximation of the temperature 
field. In the theory the extent of the potential core 
of velocity is variable, whereas in the present experi- 
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Figure 15. Axial temperature and velocity distribution 
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that the extent of the potential 
for purposes ol approximation 


ments it 
core Ol temperature 


appears 


can be re carded as constant 

Data at the conditions of | oF | 
included in figure 15, but comparison is not possible 
not investigated in refer- 
same inverse relationship is 


are also 
this condition is 
2 and 7}. The 
evident, but here there appears to be no evidence of 
a characteristic length 

2. Lateral Distributior Numerous 
traverses of temperature were made in the mixing 
zone of the chamber to investigate the lateral spread 
of temperature in the two coaxial streams. The 
generalized profiles were compared with that of the 
cosine curve as assumed by Squire and Trouncer in 
2}, and with that caleulated by Tollmien [1] for a 
free,round, incompressible jet. Tollmien also used the 
momentum transfer theory and assumed that mixing 
length was proportional to the width of the mixing 
region. However, the calculations were based on the 
assumption of a point source, and therefore in prac- 
tice similarity of profile would not be reached until 
axial distance became large compared to jet size 

Four profiles of temperature are presented in 
figure 16 to illustrate the shape over the range of 
operating conditions. In some of the runs, tempera- 
ture of the secondary air, 7: was measured just 


because 


ences 


of Ten pe rature 
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before the point of mixing, and, as expe: 
somewhat higher than the inlet temp: 
because of heat transfer from the prin 
Two profiles of temperature in figure 1 
culated by using these measured valu 
Effect of this change to 7;, was a slight s 
data points to lesser values of the temper: 
tion and to greater values of the radius fra 
At high relative velocity of the prima 
the edge of the mixing region intersects 
the 6-in. pipe in which mixing takes pla: 
tively small values of axial distance 
quently the profiles are not complete 
mixing downstream of the point of interse: 
longer in free streams, and the validity of 
with mixing in free streams is questionab 
some information can be gained by a 
and it is presented with this limitatior 
(tan axial distance of 9 r,,,, the profile of 
ture is in good agreement with the cosine « 


evel 


Son 


at greater distances there is a tendency to 
the distribution predicted by Tollmien. TI 
mental velocity profiles of Forstal] 
were a match with the 
distances studied Ratio of velocity of seco: 
primary affect. the 
temperature profile 

Although the shape of the profile tends to 
with axial distance, the magnitude of the 
half temperature r, 
of the mixing region and is plotted versus 
distance in figure 17. Theoretical results of S 
and Trouncer for velocity are included for 
Experimental findings of Forstall and S| 
well with the theory and therefor 
included in the figure 

At the condition | | 125, 7 
slowly in the region of the potential core, 
this case extends downstream 
in the theory It is t 


f 


and 
cood cosine cul 


streams doesnot shap 


is a good indication of t] 


son 


are 


agreed 


a distance of 
compared to about 8 7 
noted that at the end of the potential cor: 
perature 7, 18 1.1 73,,, and therefore the width o! 
mixing region would be about 2.2 r,., if it we 
cut off by the 6-in. pipe (width=1.7 7 Be 
the potential core the experimental data are no 
consistent, but it may be said that r, grows appr 
mately as predicted in the theory. With res 
to the two points of largest r,,, these also ar 
inally beyond the wall of the 6-in. pipe. How 
in these instances the center line of the profil 
off the center of the pipe enough to measure tl 
radii on one side of the traverse. This asymm 
has been one of the major troubles in these exp 
ments. Unknown changes in the configuratio: 
the burner, possibly due to high prevailing temp 
tures, cause asymmetry and shifts of the center 
lines of the traverse. If the curve as drawn rep 
sents the data, then the constant of proportionali' 
C in the analysis would have to be less than 0.0007 
by about one-third to bring the theoretical results 
in agreement with the present experiments. T! 
rate of mixing on this basis is correspondingly low 
than in the theory and in the other experiments 





























ns at an axial distance of about 
nsidered additional evidence that 


the extent of the 
but Yr, converges 
this point the mix- 
of the pipe (1.7 7 
refore mixing influenced by 
In this downstream region, 7, diverges 
mately as predicted by theory but again at 
values axial distance Radius of half 
im velocity was compared to that of tem- 
in three runs and about 5 
which confirms the that tem- 
is transferred at a greater rate than velocity. 
also included in figure 17 on the radius 
maximum temperature when velocity of the 
stream was low compared to the secondary 
V;.,/Vs.5=1.1 and 1.5 Convergence in 
on of the potential core is pronounced, and 
nee begins again at an axial distance of about 
In all three cases in figure 17, divergence 
11 rs.,, and this 
the chamber 


condition \ | 
is also about 11 
of about O.S5 \t 


on extends to the 


Corl 


wal 
downstream is 


¢ 
ol 
was less by 
con lusion 


aie 


haracteristic length 

runs in which V,.,./V;, was 1.1 were made to 
effect of velocity difference between the 
and secondary streams, and data from these 
included in figure 17. Curves were 
wn through the data points in these two runs 
h the mass rate of flow differed by 65 percent 


| " 
iit 


re also 





| 
| 
| 
| 
| 





by 65 


differed 
by this 


Because the rate of flow 
4 | differed 
out of four comparisons ol 
tude of V V;, has litth 
an axial distance of 187 the 
was beyond the range of the 
in one of the runs, and consequently 
plotted broken cirel 
true value 
better than 
tance On 


percent 
also amount Three 
show that the magni 
affect on the field \t 
maximum temperature 
recorder 
the radius 
is somewhat vreatel than the 
It is expected that agreement would be 
indicated in 17 at this axial dis 
the basis that r,, is not affected, it ts 
concluded that velocity rather than 
difference, determines the mixing field This same 
effect was evident in another comparison in which 
V,;, was 0.1 V Forstall and 
reached this same conclusion 

In the region of the potential core, it may be said 
that the mixing region approximates that of a two 
dimensional region semiinfinite moving 
streams differing in both temperature and velocity 
In recent years the problem of mixing in the two- 
dimensiona! with the mixing region bounded 
on one side by a quiet body of gas, has received con- 
siderable attention presented a cal- 
culation of the profile for incompressible 
flow, and much of the experimental work has been 
the determination of a scale factor ¢ to 
match the experimental with the theoretical profile 
The magniture of ¢ may 
measure of the rate of mixing 

A comparison of the temperature profile in the 
region of the potential core with Tollmien’s profile 


temperature 


higure 


ratio velocity 


about Shapiro 


betw ecn 


case 


Tollmien [1 
\ elocit \ 


based on 


be regarded as an inverse 
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of velocity in a free jet boundary is presented in 
figure 18. <A test in which | \;,» was 0.12 was 
selected as it was the nearest approach in the present 
series of experiments to the condition of still gas on 
one side of the mixing region. The other important 
difference 1s, of course, the large difference of tem- 
perature and density between the streams in these 
experiments. It appears that a value of o=—12 
brings a fair correspondence between the data and 
the curve. Liepmann and Laufer [8] report a value 
of 12 based on their own and other measurements of 
However, the conditions 


chamber has a characteristic length tha 
regarded as the extent of the undisturbe: 
core of temperature of the primary strea 
incompressible case the extent of the cor 
is variable with velocity ratio. In the reg 
potential core a convergence of the rad 
maximum temperature is noted, while un¢ 
conditions in the incompressible case t 
divergence Beyond this region, the rad) 
maximum temperature appears to have tl 
divergence. Axial distribution of temp: 
similar to distribution of velocity, and 


velocity In a subsonic jet 
differ in several respects as mentioned, and no firm 
therefore be drawn about the 
and these data are presented 
Values of ¢ 
was larger than unity, and 


pressed with reference to the temparatu 
secondary stream, an inverse relationship 
distance is obtained. This inverse relatio 
which holds without regard to velocity ratio 
general utility in the determination of the 

ture field, especially when it is noted tha 
chamber has a characteristic length. On th 
of rate of approach of axial temperature to aver 
or “mixed”’ outlet temperature, the mixing proces 
least effective when the velocities of the Streams 
equal Comparison of the field in this with tha: 
the incompressible case leads to the conclusior 
the effect of the large difference of temperat rt 
density between the streams is to retard mixing 


conclusions can 
relative mixing rate 
merely as information 
termined when | | 

these ranged as high as 18 


were also de- 


6. Conclusion 


The problems of flow, heat release, and mixing in 
an equivalent gas turbine combustion chamber have 
been examined analytically \ large primary zone 
is advantageous, but the area is limited to about 60 
percent of the total cross-sectional area Practical 
experience has also led to the use of chambers of 
about this size relationship Method of control 
currently used, although no doubt simple, leads to 
poor mixture quality in the primary zone and there- 
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fore puts a burden on the combustion process, which 
experience shows it cannot always overcome Re- 
quirements for and effects of a possible means of 
mixture control are determined analytically. Ex- 
perimental work demonstrates that the analytical 
method gives a good des« ription of the flow, heat 
release, and mixing in the controlled combustion 
chamber Although the control described would 
complicate the structure of the chamber, improve- 
ments in combustion could be a compensation. 
Study of the mixing of primary flame gas and 
secondary air, and comparison with other studies of 
mixing reveals a general similarity. Ratio of veloc- 
ity of the streams, not the velocity difference, 
determines the field of mixing, but the shape of the 
generalized profile of lateral spread of the mixing 
region is not affected by velocity. Temperature is 
transferred at a greater rate than velocity, and this 
conclusion also appears to be general in other studies. 
Some difference is noted, in that the combustion 
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Corrosion of Galvanized Steel in Soils 


Irving A. Denison and Melvin Romanoff 


The results of n 
zine as reference 


years are reported 


Phe 


mined by weight loss and pitting is interpreted in terms of the e1 
From analysis of data on 


which the specimens were ¢ xposed 
specimens having different weights of coating 


environmental conditions are suggested 


l. Introduction 


to evaluate the effectiveness of zine coat- 
itecting iron and steel from corrosion under- 
st specimens of iron and steel coated with 
thicknesses of zinc, applied by the hot-dip 
ere included in the extensive field-exposure 
ertaken by the National Bureau of Stand- 
924. The thinner coatings did not prevent 
in the more corrosive but coatings 
nominal weight of 3 0z/ft® of exposed surface 
d the development of measurable pits on 
mens at all but one of the 47 test sites for 
period of exposure |1, 2] Although the 
resented a wide range of soil conditions, it 
sidered desirable to expose specimens with 
vier zine coatings to additional soil environ- 
dominant factor, 
acidity, ete., the 
Accordingly, speci- 
same weight of coating, 
of surface were buried at 15 addi- 
tesin 1937. A sufficient number of specimens 
nized steel and of the reference materials, 
were buried to provide for removal 
specimens of each material after each of 
wriods of exposure. This paper deals with the 
of the specimens of the three materials 

for different periods up to 13 yr 


souls, 


corrosive 
was 


which 
chloride 
cause of 
ving the 
} oz/ft 


some 
content 
Corrosion 


nominal 


Zinc, 


2. Properties of the Soils at the Test Sites 


test sites represent a wide range of environ- 
conditions with respect to both the chemical 
vsical properties of the soils (table 1 For 

the hvdrogen-ion concentration of the soils 
from pH 2.6 to 9.4, and the resistivity from 
7,800 ohm-cm. At several of the test 
s are highly oxidized to considerable depths, 

in others the permanent water table is 
the surface. The environmental conditions 
al of the test sites are illustrated in figure 1. 


3. Materials 


galvanized and bare steel specimens were 
5-em) lengths of pipe 1.5 in. (3.8 em) in 
r cut from the same lots of commercial mate- 


rackets 


sites 


ndicate the literature referer the end of this paper 


easurements of the corrosior i 
materials after exposure to different soil conditions for a maximum of 13 
magnitude and progress of corrosion of galvanized steel as deter- 


of galvanized steel and of bare steel an« 


vironmental conditions to 
corrosion of galvanized 
ts for different 


the 


minimum coating requireme! 


rials. The nominal weight of the zinc coating on the 
galvanized specimens was 3.08 oz/ft?. The impuri- 
ties in the spelter in percent were: Al, 0.008; Fe, 
0.07; and Pb, 0.78. The steel, both galvanized and 
bare, was AISI No. B1010, with the nominal com- 
position in percent follows: C, 0.13 max; Mn, 
0.30 to 0.60: P. 0.07 to 0.12: S. 0.06 max. The 
interior of the specimens was coated with heavy 
grease, and the ends were closed in order to prevent 


iis 


internal corrosion 

The zine specimens were plates having the dimen- 
sions 12 in. (30.5 em) by 2.3 in. (5.8 em) by 0.15 in. 
0.4 em The impurities determined in the zine 
in percent were: Fe, 0.009; Pb, 0.095; Cd. 0.0038 


4. Distribution of the Coating on the Gal- 
vanized Specimens 


The distribution of the coating on the galvanized 
magnetic method [3] 
one 14-1n.(35 em 
random from 
on a 6-1n 
from 


specimens was measured by a 
The made 
length of galvanized pipe selected at 
the lot of specimens before burial, and 
15 em) section cut from a second specimen 
the same lot \ total of 389 measurements of thick- 
ness were made on the 14-in. length of pipe (identi- 
fied as specimen A), and 162 measurements were 
made on the 6-in. section (ide ntified as specimen B 
The measurements of coating thickness are shown 
as frequency distribution curves in figure 2. These 
curves indicate not only that the specimens differed 
widely in coating thickness, but that the distribu- 
tion of the coating on the individual specimens was 
highls The variation did not fall 


in any particular pattern 


measurements were on 


variable as well 


5. Thickness of the Outer Zinc Coating and 
of the Alloy Layer 


The thickness of the outer zine coating and of the 
alloy layer was measured by modifications of the 
electrolytic stripping methods described by Britton 
[4] and by Anderson and Manuel [5] 

The electrolytic cell consisted of a 4-0z (120-ml) 
plastic bottle from which the bottom had been re- 


Fact 
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ie cathode was a strip of zine 1.5 in 
width formed so as to fit tightly around 
riphery of the bottle. Contact with the 
made by means of a hole 0.19 in. (5 mm 
na No. 8 rubber stopper \ larger hole 
through one end of the stopper per- 
stopper to be placed over the mouth of 
The electrolyte was a solution contain- 
NaCl, 10 g of ZnSO,, and 100 g of H,O 
tion. the bottle was inverted and the test 
| in contact with the hole in the stopper 
prevent leakage of the 
battery 
initial 


ent pressure to 
By means of a 6-v storage 
a variable voltage divider, an 
25 ma was passed through the cell with 
ce as the anode, the current being read 
ding milliammeter. As the outer zine 
as removed by electrolysis and different 
zine-iron alloy and finally steel were ex- 
changes in potential associated with each 
<luced corresponding changes in the re 
rrent 
terpretation of the current-time curves in 
oating thickness is illustrated by figure 3 
mum thickness of the zine layer is measured 
orizontal part of the curve at point A. On 
of an alloy layer locally, the counter emf 
and the current proportionally decreases 
zinc is completely removed and a uniform 
the first phase of the alloy layer is exposed 
The current then remains unchanged 
cond phase of the alloy layer is exposed and 
hrough the various phases until the under- 
is completely exposed. The minimum 
s of the alloy layer, including the inter- 
phases, 1s represented by the distance BC, 
maximum thickness by the distance AD 
ported thicknesses were computed from the 
der the curve. For the present investiga- 
thickness of the zine coating was taken as 
mum thickness, that is, to point #, and the 
ss of the alloy layer was taken as 
nted by the distance RD 
maximum thickness of the outer zine coating 
the alloy layer at several points on two 
ens are recorded in table 2. These data indi- 
a large part of the zine applied to the steel 
layer 


being 


s converted to alloy 


T) kness of the oute patina, of the alloy le 


of the total coating on aalvani 2 1 spec mens 


Ir 


6. Condition of Coating, Weight Loss, and 
Pitting of Galvanized Steel, Bare Steel, 
and Zinc 


After each exposure period, one set of specimens 
was removed, returned to the laboratory, and cleaned 
free of corrosion products by previously described 
methods [2]. Photographs of the specimens of gal- 
vanized steel, bare steel, and zinc removed from the 
different sites after exposure for 13 yr are shown in 
figure 4 

The losses in weight and the depths of the deepest 
pits on the specimens, together with the percentage 
of the area of the galvanized specimens on which 
coating remained, are recorded in table 3 for the 
different periods of exposure. Except as indicated 
each value is the average of measurements made on 
two specimens 

The data of table 3 reveal that in nearly one-half 
of the total number of soils (53. 55. 61. 62. 64. 65. 66. 
and 70) the coating remained virtually intact and 
consequently provided complete protection to the 
underlying steel for the maximum period of exposure 
However, in the highly reducing soils (51, 56, 58, and 
60) and in cinders (67) little, if any, of the coating 
remained at the conclusion of the test 

In certain soils, such as 58 and 60, failure of the 
coating during the first few years of exposure re- 
sulted in a rate of the galvanized 
specimens comparable to that of bare steel. But in 
other soils, contrary to expectations, destruction of 
the coating was not accompanied by rapid corrosion 
of the underlying steel. For example, the galva- 
nized specimens in soils 51 and 56 corroded relatively 
slowly during the greater part of the exposure period 
in spite of the fact that the coating was almost 
completely removed by corrosion during the first 
few years of exposure. Even in the highly corrosive 
cinders, soil 67, the galvanized specimens corroded 
much more slowly than the specimens of bare steel 


corrosion of 


The progress of corrosion of galvanized steel, bare 
steel, and zine in the different soils is shown by the 
weight loss and pit depth-time curves of figure 5 
The values for weight loss and pitting of the gal- 
vanized specimens on which the coating remained 
intact naturally apply to the coating and not to the 
basis metal. 

The weight loss and pit-depth 
the galvanized specimens in soils 64, 65, 66, and 70 
are of particular interest in indicating that zine 
coatings were effective in protecting steel against 


time curves for 


corrosion in soils that are highly corrosive to steel, 
For example, in spite of the fact that the steel speci- 
mens in soil 66 were perforated by corrosion after- 
exposure for a few years, the coating on the galva- 
nized specimens in this soil remained perfectly con- 
tinuous throughout the entire period of the test 
In soil 70, the specimens of both zine and steel 
were severely corroded, but the coating on the 
galvanized specimens pro\y ided complete protection 
to the underlying steel 
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The weight loss and pit depth-time curves for the | occurs during the first few years of expos 
galvanized specimens (fig. 5) show that the zinc | rate decreases considerably after this 
coatings provided good protection in all of the 10 | corrosion 
inorganic soils with the exception of the highly A better comparison may be made of t 
reducing soil 56 However, in only one of the four | of specimens after each had been expose: 
organic soils (59) could protec tion be considered | Seven of the soil sites either were commor 
adequate In two of these soils, 58 and 60, the series of tests or were located in soils o 
effectiveness of the zine coating was negligible. | type, differing only in geographical locati 
Although it is not surprising that severe corrosion of possible a comparison of the two sets und 
galvanized steel occurred in these organic soils | soil conditions. The weight losses and the 
because of their corrosiveness toward zinc, it is | penetration of these specimens, both bar 
evident from the curves in figure 5 that a high rate | vanized, exposed to the same soil con 
of corrosion of zine in a soil does not indicate neces- given in table } Because the specimens 
sarily that galvanized steel would also be subject to | 1937 were not removed after exposure fo 
corrosion were the earlier specimens, it was necessa! 
polate the weight loss and pit depth-tim: 


The data presented on the condition of the coating 1088 
these specimens hig. 5 between the pomts I 


on the specimens (table 3) indicate that the speci- 
mens removed from sites 64, 65, 66, and 70 after Penen 4. Compericen between diferent tele a 
exposures for 1] and for 13 yr contained no zine | samples exposed to the same so ’ conditions * 
whateve! Because these specimens were covered yr periods 

with zinc-iron alloy layer, it might be assumed that 

the alloy layer protected the steel pipe electro- 

chemically or otherwise, but unpublished results of 

measurements of the potentials and of the polarizing 

characteristics of these materials, supplemented by 

measurements of weight loss, indicate that the zinc- 

Iron alloy laver cannot protect steel cathodically in 

soils, nor is it more resistant to corrosion than steel. 


Additional evidence that galvanic action between 
zine and steel does not fully account for the behavior 
of the galvanized specimens is provided by the data 
on weight loss and pitting of the specimens removed 
from sites 51 and 67 (table 3 Although both the 
outer zine coating and the alloy layer were completely » Site 45 is not identical wit! 


tidal marsh, is the same 


removed from the specimens during the initial 

period of exposure, the weight losses and pitting of | 13 yr. In 5 of the soils, 3 (53), 11 

the galvanized specimens during the remaining | 49 (62). and 43 (63) there is no significant differe: 

periods of exposure were considerably less than | jpn corrosion. In two soils, 23 (70) and 29 

they were on the bare steel specimens. There are | differences, particularly in penetration, ar: 

some indications that the relatively high resistance | [yn the case of soil 23 (70) the penetration 

to corrosion of the galvanized specimens is due to | greater in the 1924 specimens, but in the ot! 

an inorganic film or coating deposited by galvanic | (58), the penetration was greater in the 1937 

action between the outer zine coating and the alloy | mens 

laver or steel. The composition and properties of As pointed out previously there is a wide 

this film, which is presumably silicious in nature, are | e9ating thickness of the specimens buried 

currently under study fig. 2) and a similar but somewhat smaller 
Reference has previously been made to an earlier | the 1924 specimens. The excellent agreement 

series of field tests in which were included galvanized results of the tests of the two sets of specimens 

specumnens having the same nominal weight of coating five soils indicate that these thickness differ 

as those in the present test, namely, 3 0z/ft* [1, 2]. | were of little significance and that the differe: 

In comparing the behavior of these two sets of | corrosion of the specimens in the other two 

specimens, Blum and Brenner [6] concluded that the | must have been due to local environmental chang 

specimens in the later test corroded as a group as | cguch as drainage, aeration, etc. over the two pe! 

much in 4 yr as did the specimens in the earlier test in | of time 

10 yr. A further examination of these data, however, 

indicates that this observation is not entirely correct. 7. Effect of Weight of Coating 

The data for the 1924 specimens removed after 10 7 

yr are reported as an average rate of loss of weight Logan and Ewing [1] studied the effect of weight 

in Ounces per square foot (or penetration) per year, | of coating on the loss in weight of galvanized sp 

whereas those for the 1937 specimens are reported | mens exposed to corrosion in six soils, represen! 

as total loss of weight (or total penetration) after 4 | a wide range of conditions, for periods of 6, 8, 

yr exposure. Comparison of the two sets of speci- | 10 yr. Their results showed that the weight 

mens on this basis is difficult because most corrosion | decreased with increase in weight of ceati 
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ng value beyond which it was constant 
of this relation and other considerations, 
m was drawn that galvanized coatings 
vuund use should have a weight of not 
oz/ft? of exposed surface. Although 
Ewing showed a rough correlation be- 
t loss and pitting, no attempt was made 
rate of pitting of the specimens directly 
eight. 
to obtain a direct relation between the 
alvanized steel specimens and weight of 
depths of the deepest pits after 10-yr 
the specimens studied by Logan and 
averaged for the same six soils and 
nst the corresponding weights of coating 
¢ the values shown in figure 6, the pit 
specimens having coating weights within 
0.8 to 0.9, 0.9 to 1.0, etc. oz/ft? were 
These average values also include data 
ens of different basis metals, namely, 
h iron, wrought iron and copper bearing 
ecific effect of the basis metal on weight 
been observed. Each point shown in 
the average of measurements made on 
» 36 specimens, usually the latter. The 
the standard deviation of the pit depths 
fferent coating thicknesses are seen to 
narkedly with increase in coating thick- 
the soil conditions represented, a zinc 
having a nominal weight of 2 oz/ft? of 
surface would prevent appreciable pitting 


TABLE 5 


of steel for at least 10 yr, assuming reasonably 
uniform distribution of the coating 

As the pit depth data shown in figure 6 apply to 
specimens of galvanized steel which were exposed 


to a variety of soil environments, an optimum weight 


of coating that might be selected on the basis of 
these data might be excessive in some specific 
environment but insufficient in others. In table 5 
are shown the depths of the deepest pits on the 
specimens exposed in the earlier series of tests |1], 
from which can be determined the minimum weight 
of coating required to protect steel for 10 yr in 
specific soil environments. Data for the specimens 
buried in 1937 are also included in the table 
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It will be observed that in the group of oxidizing, 
inorgan on the specimens 
having a weight of coating of 1.07 oz/ft* or 
but no pits developed on the specimens having a 
LS] Hence it is reasonable to 


soils, pitting occurred 


less, 


coating of oz/ft 


conclude that a nominal 2-0z coating is probably 
adequate for oxidizing soils for at least 10 vr. A 


3-0z coating would provide adequate protection in 
the group of reducing iorganic except for 
highly reducing soils high in soluble salts, for example, 
These latter soils and the group of re- 
would require a 


soils 


51 and 56 
ducing organic 


soils and ( inders 


coating in excess of 3 oz/Tt 


8. Summary 


This report contains the results of measurements 
of corrosion made on specimens of galvanized steel, 
bare steel and zine exposed to different soil conditions 
Che nominal weight of the 


for periods up to 13 vr 
but the actual thickness 


zim was 5 oziit 
varied over a wide range 
of the 15 soils to which the specimens were exposed, 
the coating on the galvanized specimens remained 
virtually intact throughout the entire duration of the 
field test. The galvanized specimens were especially 
resistant to corrosion in alkaline soils that were highly 
corrosive to bare steel In only two both 
organic, failure of the zinc coating after relatively 
short exposures was accompanied by marked corro- 


of the The high 


coating 


soils, 


sion steel corrosion 


shown by most of the galvanized specimens after the 


In approximately one-half 


resistance 


outer zine coating had been virtually 
corrosion is tentatively attributed to a { 
ing, probably silicious in nature, that wa 
cathodically by galvanic action betwee) 
zine coating and the allov laver or the 
steel 

The results of these tests show that th 
weight of zine coating required to protect 
corrosion for a minimum of 10 yr dep¢ 
nature of the soil environment. In thes 
tests, a 2-0z coating was sufficient to prot 
inorganic oxidizing soils, but for inorgani 
reducing soils a 3-0z coating was require 
reducing, inorganic and organic soils requi 
in excess of 3 oz/ft 


The field tests described were initiate: 
stalled, and until 1946 were conducted 


supervision of K. H. Logan 
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Ficure 1 Environmental conditions at typical test sites. 
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ries clay at El Vist t site 70. Merced silt k at Button 
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Distribution of coating on test specimens A q , coating ¢ 


4; length, 14 in.; number of observations, 389. @ Specimen B 
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Corrosion of Low-Alloy Irons and Steels in Soils 


Irving A. Denison and Melvin Romanoff 


The 

14 soils for 

as determined by 

of the materials and the 
were exposed. 


results of measurements of the 


per ods ip to 
we 
of the 


nature envire 


l. Introduction 


National 
with number 
pipe materials 


Bureau of Standards 
of manufacturers and 
initiated an extensive 

of the of certain materials 
derground construction. Examination 
specimens removed during the early 
«posure indicated that certain alloy steels 
ed resistance over ordinary 
teels. Consequently, specimens repre- 
reater variety of alloy irons and steels 
ed at at intervals during the 
he exposure tests. The condition of the 
removed from the various after 
longer periods of exposure has been 
of reports, and in 
report on the specimens buried in 1932 
hed [1] The present report contains 
of measurements of weight loss and pitting 
of 10 of 


he 


in 
a 


corrosion 


corrosion 


the sites 
sites 
a series 


progress 


varieties irons and steels that 


sed in 1937 and removed from the test 


five periods of exposure, ranging from 
duration. The properties of the 
sites, the installation of the specimens 
and the methods in cleaning the 
removed prior 1946 are described in 
vers [1, 2]. Specimens subsequently re- 
re cleaned by immersion in fused sodium 


soils 


used 


Lo 


Description of the Materials 


mens were in the form of plates 12 in 
wide, ranging in thickness from 


25 in. 


‘om posi 


13 years are giv 
loss and pit-depth measurements are correlated with the composition 


; 


corrosio! 


mental conditions to 


on 


of 10 low-alloy 
rhe magnitude 


irons and steels exposed 
and progress of corrosior 


which the test specimens 


0.175 to 0.265 in 
are given in table 1 


The compositions of the materials 


3. Results of the Exposure Tests 


The corrosion of the different materials in typical 
soils is illustrated in figures 1 and 2. Each of the 
four environmental groups based on aeration is rep- 
resented as follows: good aeration, soil 55; fair aera- 
tion, soil 66; poor aeration, soil 61; and very 
aeration, soil 56 

The extent of corrosion damage was determined 
from the loss in weight of the specimens as a result 
of the exposure, and from the depths of the deepest 
pits. The results of measurements for the 
different periods of exposure are given in tables 2 
and 3. Unless indicated otherwise, these values are 
the of measurements on two specimens 
The exposure periods did not differ from the average 
values given in the tables by more than 5 percent. 

A number of the specimens exposed for the longer 
periods in the more corrosive soils were perforated by 
pitting. Because the depths of the pits producing 
these perforations would obviously have been greater 
after successively longer periods of exposure if 
thicker specimens had been used, some adjustment of 
the values for pitting of the perforated specimens was 
desirable. An approximate adjustment was made by 
multiplying the average penetration of the perforated 
specimens by the pitting factor, defined as the ratio 
of the maximum to the average penetration. Deni- 
son and Hobbs [4] showed that the pitting factor of 
plain irons and steels decreased with time during 
the early periods of exposure underground but 
became approximately constant after exposures for 6 
8 vr. The value of the pitting factor was con- 
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soluble 
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serated heavy clay containing a high concentrati 
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FiGuRe 2 Corrosion of ¢ 
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76 Cr, 0.43 Cr. Soil 55, well oxidized 


u-P steel, 1.02 Cr, 0.42 Cu 
! ient in soluble salts 
t of soluble material 


serated loam containing a high concentrath 
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sidered by Denison and Hobbs to be characteristic The average depths of the deepest p 
of soils and Logan [2] showed that it varied directly | materials, given by the values of the 
with the aeration. Adjusted values of pit depths y-intercept), indicate that the low-alloy 
recorded in table 3 are enclosed in parentheses greater initial pitting rate than the 
The over-all behavior of materials and the effect of However, as the exposure increased, 
the various alloving constituents on the corrosion of pitting of the alloy steels diminished m 
iron and steel in the soils is indicated by the weight than the rate for the plain steel so that 
loss and pit depth-time curves shown in figures 3 and | the order of the materials was reversed 
{ In preparing these curves, the values for weight The weight loss and pit depth data for 1 
loss and pit depth for each material in all of the soils, | chromium and chromium-molybdenum 
except 51,? were averaged for each period of exposure KK, D, E, and H, (fig. 4 and table 4) exh 
The logarithms of these average values were’ then | but somewhat greater effects of the al 
plotted against the logarithms of the periods of stituents than the copper-molybdenum a 
exposure nickel irons and steels. Chromium rei 
Che curves shown in figure 3 and 4 conform to the weight losses in a fairly regular manner, but 
equations the initial pitting rates of the steels Ho 
rates of pitting of the allov steels deer 
rapidly with time than the rate for plain 
The separate effects of chromium and mo 
on the pitting of steel are difficult to « 
because the chromium steels also genera 
molvbdenum It will be observed that st 
D, containing | and 5 percent of chromiun 
T+loe k 2 tively, had pits of about the same depth. W 
this range chromium alone in excess of 1 or 2} 
loc W=—4x loo low k’ does not appreciably increase the resistan 
n “ . material to pitting. On the other hand, th 


Hence n(u) is the slope of the line and k(k’) is the of molybdenum in reducing pitting is quit 
because all of the chromium steels cont 


intercept on the log P(W) axis =, ~ 
Equations 1 and 2 were derived originally by molybdenum, KK, E, and H, had shallower pits 

Logan, Ewing, and Denison [5], and by Martin [6], the straight chromium steel D 

respectively The pit-depth—time curves, figures 3 and 4 

cate that the change in the rate of pitting wit! 

initial corrosion rate of the materials and the change depends on the composition of the steel. In ger 


in the rate with time, were calculated by the method the greater the amount of nickel, chromium 
By molybdenum in the steels, the greater was the cl 


in the rate of pitting of the steels with time 

be assumed that the alloving constituents 
the formation of corrosion products, which 
effective in diminishing the rates of pitt 
alloy steels with time. The results of a 
study of the effect of corrosion products on ( 
of pitting of steel [8] suggests a probable expla 
for the observed differences in corrosion rats 
that study it was observed that in soils in w 
corrosion products of ferrous metals diffused 
into the soil, the rate of pitting changed | 
little with time, but in soils in which the « 
products remained in close contact with the c 
surface, the rate of pitting diminished w 
often becoming negligible after relative 
periods of exposure 


where P is the depth of the ceepest pit at the time 7 
and W is the weight loss at the time 7 
Converting to logarithms 


The constants of the equations, expressing the 


of least squares for each material in each soil 
means of these constants, values of the average 
weight loss and pitting of eac h material in all of the 
soils were calculated for the maximum periods of 
exposure. These values, together with the constants 
of the equations and their standard errors, are tabu- 
lated in table 4 

In order to estimate the probability that the weight 
loss or pitting of each material was significantly 
different from the corresponding values for plain 
steel, taken as the reference material, the standard 
t-test was applied, and from the calculated values of 
t, the probability of the differences being due to 
chance was obtained 7] 

The weight loss and pitting data for the coppel 
molybdenum open-hearth irons O and N and for the 
nickel-copper steels J and B, presented in figure 3 . . 
and in table 4 show that these allovs corroded slightly | The conclusions that can be drawn from 
less than the plain steel A. However, the probability | in table 4 and in figures 3, 4, and 5 are ne 
is high that several of the observed differences are limited because these data indicate only the 
due to chance. Because mill scale was not removed | behavior of the materials under a wide Vv 
from the specimens of steel B before burial, the data | e@Vironmental conditions. Whether a p 
recorded for this steel cannot be taken as truly | alloy is more resistant to corrosion than plain 


representing the effect of additions of 2 percent of im a specific soil environment obviously cal 
nickel and 1 percent of copper on the corrosion of predicted from these average values. Altho 
steel in these soils average rates of pitting of the alloy steels « 
git arora more rapidly than the pitting rate of plain 

«0 periods only ere — ee | should be recalled that this effect is a cons 


3 The data for the specimens 
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mment as well as of the material For 
a very poorly aerated 
ild tend to diffuse and migrate outward 
rather than to form protective lavers 

conditions, the higher 


s Under 
| pitting of the illoy steels would tend 


soil. corrosion 


such 


ined 
the materials 
the 


according to 


yarison ol the behavior of 
ent environmental 

for the classified 
bles 2 and 3 calculated on a 
sis for each all periods ol 
he weight losses and pit depths of the 
eel A being taken as 100 percent The 
these values for all periods and for all 
same environmental groups are given in 
cause the data for soil 51 are incomplete, 


not included in cal itlating the relative 


conditions, cor- 
soils 


were 


material for 


ve corrodibilitv of the 


is measured by weight loss was not affected 

bv differences in the aeration of the soils 
the 4- to 6-percent chromium steels wer 
ted in all of the very poorly aerated soils, 

rs. These steels had higher initial rates 
than plain steel, but the rates decreased 
vy with time in most of the soils. In the 
ated soils, however, this high rate of pitting 
throughout the exposure period, probably 
ynditions were not favorable to the forma- 
which would have diminished the 
is the behavior of these 
nders (table 3 and fig. 5), where the initial 
all the steels was about the same The 
the plain steel continued at a high rate for 
exposure period, but most of the pitting of 


low alloy 


irons 


bercles, 


te In contrast 


§-percent chromium steels occurred during 
t vr of the test, there being only slight 
n pit depth of these steels after that time 
ven more marked in the case of the steels 
x molybdenum in addition to chromium i 
n which there was practically no additional 
ter the first 4 y1 
containing coppe! molybdenum also 

resistance to pitting in cinders after 4 yr 
ire than the plain steel, although to a lesser 
in the 4- to 6-percent-chromium steels 


and 


4. Summary 


measurements 
irons and 


contains the results of 
several low-alloy 


port 
on made on 
exposure to different soils for periods up 
Empirical equations fitted to weight loss 
pth—time curves permitted the initial rates 
loss and pitting and the change in these 
h time to be evaluated The general effect 
ving elements was to reduce the initial rate 
mn of the alloys as measured by weight loss 


but to increase the initial rate of pitting Except in 
very poorly aerated soils, the rate of pitting of the 
alloy steels diminished more rapidly with time than 
the rate of pitting of plain steel, with the result that 
the maximum depths of pits after the maximum 
period of exposure were less on the alloy steels than 
on the plain steels 
the most pronounced effect on weight loss 
maximum reduction in pitting, molybdenum also was 

Chromium and molybdenum were pat 
effective in reducing the 
steels exposed to cinders 


(Chromium was observed to have 
Dut tor 


hecessary 


ticularly 


Corrosion ot alloy 


Che field tests described in this paper were planne al 
and installed, and until 1946 were conducted unde 
the supervision of K. H. Logan 

The writers gratefully 
tion of O. B. Ellis, Research Laboratory 
Co. and of C. P. Larrabee, Research 
United States Steel Co. in making available 
equipment for cleaning the specimens by the 
hydride process 
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Loading of Quartz Oscillator Plates 


Leland T. Sogn and 


the surface 


investigated 


loading 
experimentally 


plates have bee 
applied to different 
fited only 
FF 
the metal was applied 
In both 
Activity-frequency 


roved 


when the 
g an extensional mode 
near the 


ut crvstals employvir 


types of 


characteristics Over a 


ip 


l. Introduction 


tric quartz crystal units are extensively 
frequency controlling and _ filtering 
electronic circuits of their ex- 
onstant equivalent electrical parameters 

However, coupled with the fundamental 
there are secondary modes of 
that, under conditions favorable to them, 
ome active that thev affect both the 
nd frequency of the fundamental control- 
as well as other critical parameters of the 
iivalent circuit Various methods, such 
ntouring and edge dimensioning, have been 
ninimize the effects of this coupling. Ex- 
: recently concluded at the Bureau indicate 
tive physical loading on the piezoelectric 
te surface gives approximately equivalent 


as 


because 


mode 


su) 


2. Preliminary Experiments 


ations made on 0.5- by 0.6-in. hand-finished 
ates a few years ago vielded many interest- 

that pxpvided the basis for the present 
tion. When the electrodes coupling the 
to the oscillator circuit were placed on the 
portion of a wedge-shaped plate, not only 
frequency lower but the activity was higher 
en they were placed on the thinner portion 
ate. In order to investigate this phenome- 
: eflectively, a small brass electrode \ in 
ter at the base was constructed. This elec- 
ule it possible to couple much smaller por- 
he quartz plate to the oscillator circuit and 
the effect of coupling to local irregularities 
ystal surface that resulted from the applica- 
onuniform finger pressure during grinding. 
ted, it was possible to excite a large number 
neies by coupling areas of different thick- 
oscillator circuit. The previously men- 
nverse frequency-activity relationship was 
t and was found to be very consistent as the 
was moved about on the surface of the 
was noted that whenever a certain fre- 


he 


T+ 


parts of the crystal surfaces 
metal was applied on or near the ac 


periphery of the 
oscillator plates loading increased the 
temperature range 


Philip A. Simpson 


illator 
were 


of two types of piezoelectric 
Various amounts of 


Thi 


juart 

Wood's metal 
kness-shear-mode plates were bene 
central area of the 
responded favorably 


but 
when 


tive plate 
of vibration 
plate 

crystal unit veral 


HO° tt 90°C were 


fold 


greatly 


quency was excited by coupling to different areas of 
the crystal plate the activity was the same. It was 
also noted that the lowest frequency was associated 
with the highest activity 

Similar studies were made later on machine-lapped, 
slightly convex quartz plates. Results were much 
the same, except that random distribution of fre- 
quencies was replaced by orderly progression as the 
electrode was moved from the center along a diagonal 
toward a corner of the plate. Areas to which the 
oscillator circuit could be coupled in order to excite 
a certain frequency were arranged in the form of con- 
centric rings about the center, with the lowest fre- 
quencies near the center and the highest at the 
corners 

The frequency associated with the highest activity 
was not the lowest frequency obtained by coupling 
to the exact center of the plate, but was a slightly 
higher one excited by coupling to an area immediately 
surrounding the center. In a plate whose frequency 
increased 25 ke between center and corner, the fre- 
quency with the highest activity was 4 or 5 ke higher 
than that of the center. It was this fact that sug- 
gested the hypothesis that for each crystal plate 
there is one frequency and electrode position that will 
give the highest activity, and that physical loading 
of the plate would be beneficial through causing it to 
oscillate nearer this optimum frequency. 


3. Detailed Experiments 
3.1. Concave Crystals 


The first crystals loaded were a group of 10-Me 
plates, 0.6 by 0.5 in., with their centers approxi- 
mately 0.00005 in. thinner than the corners. The 
activities (rectified oscillator grid current) of these 
plates unloaded varied between 0.0 and 0.2 ma over 
the range—60° to +90°C. The frequency curves 
also were very poor, showing many irregularities 

The plates were then etched 3 to 5 ke and loaded 
back the same amount so as to maintain the same 
relationship between their fundamental and secon- 
dary modes. The load was Wood’s metal arranged 
in a ring concentric with the center of the plate and 
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approximately ‘4 in. in diametet It was applied 
with several light strokes rather than one heavy one 

The activity of these plate s increased to 0.7 ma o1 
revealed 


tests 


the 


and subse quent temperature 


ire quency 


more 
that the 
were much smoother 

In order to investigate further the effect of loading 
several 10—M« etched loaded in 
approximatel The showed 
immediate improvement continued to improve 
until lowered 6 ke by loading, after which further 
loading resulted in a slight loss in activity (fig. 1 
\Nloreover 
dition 
the temperature range; the fl 


hg. 2 


curves, as well as activity 


eurves 
and 


crvstals were 


1—ke steps ervstals 


which in the unloaded con- 
throughout 


the ervstals 
worked intermittently, oscillated 
tuations decreased in 


number and in siz 


3.2. Flat Crystals 


AT-cut plates, 0.6 by 0.5 in., of a lower frequency 
also were prepared and showed a favorable response 
to loading The plates whose frequene les were 
,.700 and 5,500 ke, had a very slight convex 
contoul less than 0.00001 in They were etched in 
5-ke steps and loaded back to the original frequen \ 
The increase at and best re- 
reducing the frequency 
loading resulted in a 
the had the 
it created an area 


around 


activity began to once 


sults were attaune d with a load 


Further 
Evidently 
Sane that is 
of low frequency in the center of the plate 

Several different loading were 
results were obtained when the major length of 
the load was parallel to the X-direction, although a 
circular pattern did give good results with concave 
crystals. Besides the circular pattern, an elliptical 
load, wide-band load, narrow-band load, and pinhead 
load were used (fig. 4 \ narrow band applied 
paralle! to the Z’-direction appre ared to have a dele- 
terious effect on the performance. The most favor- 
able pattern appeared to be the elliptical one, al- 
though the wide-band load also improved the crystal 
to a large degree. The narrow-band load in the X- 
direction produced only a slight improvement, which 
neve approached that of the elliptical or wide band 
In general, it appeared that a load was most beneficial 
when its shape conformed to that of the active central 
area of the oscillator plate 

Interesting results were obtained with a pinhead 
size load. Initial application of load decreased the 
frequency as expected. However, after the crystal 
fundamental had been reduced about 0.2 percent, 
the frequency suddenly jumped to many kilocycles 
above the fundamental of the unloaded crystal. Ex- 
plorations with a small probe electrode revealed that 
the plate was inert along a line through the center 
parallel to the X-axis, and that each half of the plate 
was oscillating. Lycopodium powder patterns fur- 
ther confirmed this. 

Measurements were made on 0.6- by 0.5-in. AT- 
cut plates with the X-direction parallel to the long 
dimension. Their frequency was approximately 


5 ke hg ) 
Activity 
effect as contouring 


decrease in load 


patterns of used 


Jest 


6,702 ke. Elliptical loads differing in 
width, were used. The plates were 
then loaded back to their original frequer 


not 


maintain as closely as possible the origi 
ship between primary and secondary mov 

Associated with the general improv 
crystal unit there is, as expected, an ine: 
(), or quality factor. A crystal unit ma 
sented as a tuned circuit as shown in figw 
C, is the equivalent series (motional) ¢; 
the series resistance, 1, the equivalent se1 
ance, and C, the total effective capacity sli 
resonant circuit. The Q of this « 
then be determined by the formula 


series 


i) 
, tr A fhe 


where Af j j the difference betwee 
resonance and resonance of the crystal un 

The loading of the crystal plate affect: 
R, very much but produced, at most, a 
change in C Both Af and R, tended to d 
first, resulting in a rise in Q. At 
in the @ occurred, probably due to coupling 
the fundamental and some other mode asso 
with the oscillator plate. However, th 
unit ( was still several times the original o1 

When the frequency change produced by 
reached approximately 1 percent of the fundam 
R, began to increase, resulting in a drop in 
However, the Q continued upward because A/ 
still decreasing (fig. 6 Further loading cause 
frequency to jump to a value several kilocyecles 
that of the unloaded plate (fig. 7 It appears 
the loaded area became so far off frequen 
thinner region farther from the center becan 
In other words, this was just another cas 
previously mentioned “‘pinhead jump.” 

The area of the load also had a definite « 
the Q. For a load one-third the length 
X-dimension and 0.3 the width of the Z’ 
the 6,702-ke plate rose steeply to a high va 
remained high until shifted 65 ke. Attempts 
shift 70 ke resulted in a frequency jump. For 
longer in the X-direction, the Q rose less st 
but heavier loads could be applied befor 
frequency jumped. In all the (’s 
approximately the same maximum value. A 
crease in the width of a given load was a 
panied by an increase in Af, but its effect 
was unpredictable. 


some 


cases, 


3.3. Off-Center Load 


In order to determine the effect of placing th: 
off-centet, 12 BT-cut plates were prepared at ap| 
imately 4,930 ke. These crystals were 0.6 by 
in., with the X-axis parallel to the short dir 
and had a slight convexity. Areas \, \¢, 0! 
diameter were then plated on both sides w 
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ocated as shown in figure 8 Contact was 
with the plating so that the 
as well as a In order to map 


a pencil point under moderate 


gold served 
rode load 
ve area 
is then moved from the edges toward the 
| all The 
iow the position of the pencil when this 


The 


oscillations ceased points in 


region W 1 was considered the 


finite effects were noted Although the 
were not over the center of the erysta! and 
was therefore apparently outside the field 
the central portion of the 


center ot the 


\ the electrode 

still approximately in the 
a. However, the 
also encompassed 
le and usually extended bevond it 
on of the load did shift the region of 
n the very much, but did 
he active area when plac ed off-center 


active area was shifted 
some of the region unde 
fig. S 
not 


central area 


3.4. X-cuts 


oscillator plates have their 
acteristics if their contour is concave rather 
This type, which vibrates in an ex- 

compressional is more difficult to 
o oscillation if the contour is flat or convex 
he resulting crystal units have low Q 


etore 


best pel for mMm- 
vex 


mode 


seemed logical that to obtain a rood 
v loading, the metal should be applied neat 
; of the face rather than in the center, and 
ts justify this reasoning 


} 


Hates were approximat ly in. square and 


mtour such that the corners were 0.0002 in 
than the center They 
le air-gap holder, and the air gap was ad- 
or MaXimum activity 
ise of the contour all of the plates had fairly 
aracteristics before loading One unit X 2, 
inloaded Q of approximately 100,000, but 
nereased to approximately 262 OOO bv the 
load The loaded area band 
de along the edges, and the material was 
metal (fig. 9 Table 1 shows the resulting 
ement in this unit. The small changes in (¢ 
result of varying the air gap to keep the 
maximum 


were mounted in an 


is of v was a 


ata 


area by w dening the load had 
area of the load in 
plates It re- 
had no predictable 


Slightly ine reasing the 
the same effect as increasing the 
the 
sulted Im an 
effect on R 


thickness-shear-mode 
increase in Af but 


Cast of 


4. Summary 


The 
and condition of the oscillator plate but also with the 
Size, load The first effect 
observed in heavily loaded flat thickness-she aul mode 
plates, due to the establishment of a 
large frequency differential between the 
center and edges, causing the lower frequency vibra- 


effects of loading varv not onlv with the tvpe 


position and mass of the 
is probably 
resonance 


tions originating in the central regions to be almost 
completely attenuated hig 10 Bec ause 
ondary fundamental! are con- 
fined to the e« ntral region and are ¢ ffectively reduced 
in amplitude, the peripheral area of the plate is 
unable to and remains relatively inert 
This concentration of activity in the center isolates 
the vibrating quartz and the and prevents 
reflections ol energy from the crvstal edges 

The effect of loading observable in coneave thick 
ness-shear-mode plates and X cut plates operating 
extensional mod parallel to the thickness 
one primarily of dimensioning 
The 
purpose of slightly changing the frequency of certain 
plate so that the distribution of like 
frequency areas over the 
favorable for oscillation 


the Ser 
modes coupled to the 
oscillate 


mount 


In an 
appears 
certain areas by 


also to be 
loading load here serves the 
regions of the 
plate as a whole is more 
In concave plates this is 
accomplished by the addition of the load to the 
central regions. It that thickness-shear 
mode erystals of this contour vibrate at a lrequen \ 
the optimum Irequency of the plate as a 
because the central controlling 
The loading of the central area makes its 
frequency slightly thereby establishing con- 
ditions more favorable for the plate as a whole 
Flat or convex X-cut plates vibrating in extensional! 
mode in the thickness direction are improved by the 


appears 


above 
whole 
thinnest 


area is the 


lowe I 


application of the load to the peripheral areas of the 


X faces 
surrounding quartz having a higher natural frequency 
than that of the central areas of the plate Loading 
the outer areas tends to correct this condition by 
reducing the frequency of the higher frequency 
peripheral areas. This loading may be a complished 
by applying a foreign material to the outer portions 
of the plate, as illustrated in figure 9, or by hollow 
grinding the edges, as shown in figure 11. In the 
latter case, the overhanging quartz may be thought 
of as constituting the load. This type of loading 
effective with plates of relatively low 


This improvement is possibly due to the 


Is more 
frequency 
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In provement of a concave HT-cut oscillator plate CC-8 as a result of loading. 


oaded run is due to the melting of the Wood's metal, 
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Effects of Moderate Biaxial Stretch-Forming on Tensile 
ind Crazing Properties of Acrylic Plastic Glazing 


B. M. Axilrod, M. A. Sherman, V. Cohen, and I. Wolock 


1. Introduction 


gh polymethyl methacrylate glazing in air- 
equently prepared by a forming process that 
s the material, there is little information 
on the effect of this stretching on the tensile 
material. Some data 
Northrop Aircraft 


specimens 


ng properties oO! the 
tvpe were obtained at 
Tensile 
om pleces of polymethyl] methacrylate that 
n stretched uniaxially The pieces were 
ed about 60 percent while at 129° C (265° F 
oled while held at this elongation. It was 
hat at both room and subzero temperatures 
cimens oriented transversely to the direction 
ch were appreciably weaker than the longi- 
specimens ; temperature the 
specimens showed appreciable permanent set 
ast to the former 
effects of hot stretch-forming on polysty rene, 
rial somewhat similar to polymethyl meth- 
in forming behavior, have been reported by 
|. It was found that uniaxial stretching of 
hundred percent greatly increased the tensile 
th, the elongation at failure, and the “crazing 
h in the direction of stretch; the tensile 
gth was greatly reduced perpendicular to the 
on of stretch. Also the tensile strength of 
hot-stretched first longitudinally and then 
ersely roughly 200 percent was reported as 
increased for both directions. 
experiments described in this report 
to gain more information on the effect of form- 
the crazing and other properties of polymethyl 


tests were made on 


also, at room 


} 


were 


vestigation was conducted under the sponsorship 

mimittee for Aerenautics; tt report is a condensation ¢ 
CT N2779. 

brackets indicate the literature references at the end of t 

¢ [2] the test co cated nor was i 


razing strength”’ was a str crazing t 


of the National 
f NAC 


nditior were 
g or a solvent 


determined on both 
of sheet material 
elongation, strain 
threshold 


resistance to 


The 
unformed 
strength 
the onset of 
crazing 


methacry late 
formed 

include d 
and 


stress tor 


propertl = 
and 
tensile 


pieces 

total 
CTaZing 
and 


stress al stress 
stress-solv« nt 
weathering 

The major portion of the work was done on mate- 
rial stretched biaxially to an elongation of about 50 
that is, about 50 percent in all directions in 


percent 
A few experiments were also 


the plane of the sheet 
made on pireces stretched shehtly about 7 to 20 per- 
This work was carried out as one phase of a 
concerned with factors affecting 


properties ol laminated 


cent 
researe h 
the crazing 
acrvlic glazing 


program 


and strength 


2. Materials 


commercial cast poly- 
in thick- 


The materials used were 
methyl methacrylate sheets 0.12 to 0.15 in 
ness. The samples used for all experiments except 
the exploratory work were obtained directly from the 
manufacturers and were masked on one side only, as 
to make laminated acrylic 


which included both the 
‘ 


is done for sheets used 

glazing. These samples, 
general-purpose grade and the heat-resistant grade, 
18- by 0.15-in. sheet from 
They are referred to 
samples and are 


consisted of one 36- by 
each of three production runs 
subsequently as “representative 
identified as follows 





3. Apparatus and Procedure 


3.1. Forming Process 


a. Equipment and Procedure 


is that would produce 
10 in. in diameter 
according to sugyvestions offered by 
of the Rohm & Haas Co 

diagram of the forming equipment ts 
sheet of 
State Is 
forming 
vessel 


\ vacuum 


flat biaxially stretched disks about 


lorming apparat 


wis designed 
W. F. Barto 

\ schemati 
shown in higure | In this 
material, A, heated to the 
flange of the evlindrical 
is created in the 


apparatus a 
acrvli rubbery 
clamped to the 
vessel, B. A partial vacuum 
by connecting the latter to an 
The differential is controlled by 
valve, ( unclamped part of the 
The form D, an open-ended 
eviindrical tube a little smaller in diameter than the 
forming vessel and constrained by the ruide E, is 
inserted into the vessel The pressure differential is 
then removed quickly by admitting air through the 
plug valve, F, so that the stretched acrylic sheet 
shrinks about the end of the The sheet cannot 
retract completely; th portion remains 
uniformly stretched across the open end of the form 
The formed aervlic shaped like a top hat, is 
temperature before 


evacuated tank 
the plug 


pressure 
sheet is 


and the 
drawn into the vessel 


form 


central 


she et 


cooled to room in the vesse] 


removal 




















TO VACUUM 
—— TANK 


Se he matte ium-forming appa atua 


sheet to be formed forming vessel: ( lve to evacuated tank 


guide; F 


valve t 


In practice, the forming operation is 
quickly as possible so that the acrylic sheet 
be in the rubbery state when the pressure d 
the time removal of 
from the oven until forming is complete is 
1 min. The forming apparatus, 
drawing of a sheet and disassembled for thy 
shown in figures 


Is removed: irom 


assemble« 


drawn sheet is 


respective ly 


¢ 
O1 il 


b. Uniformity of Forming and Equation for Elong 


the amount of 
face of the biaxially 
pieces following experiments 
10-in Lucite HC201, which 
axially stretched 150 percent after heating 1 
was marked off in 1-in then heated t 
and allowed to assume its original size | 


whether 
the 


To determine 
is uniform over 
the 


disk of 


were 
had 


squares 


on the resulting disk were still equidistar 
the standard deviation of the 
ment and marking errors, indicating that the 
of stretching was reasonably uniform ove 
Next, another piece of Lucit 
marked with a grid having a 
of 15 mm; the piece was then biaxially hot-st 
to an elongation of 150 percent. The lines 
flat top of the stretched dome were still equ 
to within +5 percent, verifying that the st 
was reasonably uniform 

The formula used for calculating the am 
a formed disk is 


) percent 


of the disk 


was square 


biaxial stretching in 


Figure 2. Vacuum-forming apparatus with a sheet of 


plastic in place ready to be formed 
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, } 
} 


ackground 


rraphic record I 


e=100(yt,/t;—1), l 


are 
This 
of the 
constant on stretching 
verified by 
piece of both 160- 


the elongation in percent, and ¢, and ¢ 
and final thicknesses, respectively 
; based on the fact that the volums 
emains essentially 
operty of the 
the density 
med 
samples 


terials was 
1] 


it 


na 
of a sma 
and unformed material from each of 
The of Lucite HC201 
lecrease in density of 0.8 percent as a result 

to 160-percent The other thre« 
howed density changes of less than 0.2 per- 
T 


CSsuUul 


sample 
strain 


of this amount of stretceh-forming 


3.2. Standard Tensile Test 


mad in most 
1011 of Federal 
tests the threshold 
by 
ediately applied a sudden momentary pres- 

sensitive cross head of the testing machine 
. jog in the load-extension record drawn by 
The strain and the stress 
reshold could thus be readily obtained from 
1. The observer viewed the crazing against 
using north davlight or fluores- 


tensile test wert 
accordance with Method 
on L—P—406a In these 
razing was noted 


il dard 


visually an observer 


it 


2 


Vacuum-forming apparatus partly disassembled 


after form ng a sheet 


the forn 


p aati 


ece is on the end o 
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3.3. Stress-Solvent Crazing Test 
In the stress-solvent crazing tests, tapered tensil 
specimens of the same dimensions as the long-time 
tensile test specimens described in section 3.4 
placed undet load ih a hydraulic machine 
benzene applied, and the load maintained lor 4 min 
From preliminary trials on other specimens, the loads 


We 


rt 


testing 


on these specimens were selected to produce crazin 
over a part of the tapered portion The 
mens of each formed or control piece wet 
shghtly different in an effort t 
threshold at different parts ol he t 
the specimen 

Benzene was applied to the 
portion of t 


two speci 
e teste d W itl 
the 


apere ad portion ol 


loads Nn 8) locat 


DY il 
Subse que ntly 


unc I 


central 
he specimen with a brush 
the solvent crazed 


suitable | 


were eXamuined 
t of the « 


specimens 


ivhting and the exten noted 


razing 
3.4. Long-Time Tensile Test 


The long-time tensile loading cabinet for test 
high relative humidity about percent 
shown in figures 4 and 5. A similar cabinet without 
a front cover or blower and with an interior instead of 
exterior light was used for tests at 50-percent relative 
humidity, the condition in the controlled atmosphe rm 
in which the cabinets 

In each cabinet four specimens can be tested simul 
taneously. The load on the specimen is applied by 
a 300-lb-capacity weigh beam through a turnbuckl 
A pair of alinement holes in each end of the specimen 
and an alinement hole and an alinement pin in the 
clamp, as shown on the specimen at the left of figure 
4, facilitate alining the clamps and specimen. The 
specimen alinement must be done very carefully, as 


ing at 
Q5 


room were located 


otherwise stress crazing occurs much sooner on one 
face than on the other or along on« 
across the width of the specimen 

The humidity cabinet (fig. 4 that 
directs air against cloth wicks dipping into a tray of 
water. The relative humidity is readily measured 


edge rather than 


has a blows I 


Interior of 


lfesting vf gh ré 


FicuReE 4 long-time tensile loading cabinet used fo 


ative humidity 





L.-P-406a, except that it was carried out f 
instead of 240 hr, the recommended time 

Light-transmission and haze measureme 
taken before and after the weathering test 
pivotable-sphere hazemeter, following the 
in ASTM Method D 1003-49T 

To permit the measurement of shrinkags 
were ruled on each specimen about 2 tn. apa 
the weathering The scratches were measu 
a steel rule vraduated to hundredths of an 
measuring the distance between scratches 
power magnifier was used, the distance be 
mated to a thousandth of an inch 


3.6. Degree of Forming of Representative 


Exploratory tests were made with pieces o 
plastic sheet biaxially stretched slightly, 7 
percent elongation, and moderately, about 
cent. The latter elongation is an amount 
be attained at some locations in formed 
, , enclosures [3,4] The results of tensile an 
with wet- and dry-bulb thermometers placed neat solvent crazing tests indicated that the craziz 
the exhaust part of the blower The relative humid- erties peta pA threshold of ahaa crazing 
ity Is maintained at 95 2 percent standard tensile test and threshold stress rm 
To avoid heating of the cabinet, the fluorescent solvent crazing, were unaflected or only very 
lights used for observing crazing are placed just out- aftected - biesiel atestching T te 30 ‘ 
side a window in the top of the box. Mirrors However for the 15-percent-stretched _ pix 
mounted behind ead slightly above the specunens cTazing properties were considerably chan ed 
direct the light against the latter, which are viewed eon age valk oe pereamcgge seemanan 
against a black-felt background crazing appeared to increase about 50 
Che tensile specimen used for the long-time loading A oe ordinals it was decided to form pieces fre 
oma a a -long reduced —e tapering UNI- | choot of the representative samples to an elonga 
formiv in width from 0.50 tn. at the maximum cross , ~ eat ? 
section to 0.33 in. at the minimum In this way the “ cogs mcg. Age! sooner each sheet 
stress in the reduced section decreases from a value eng 7“ f ' sed was heated in om oven to 
of S, at the minimum section to 2/3 S, at the other onto Hood of 1an° es 140° C. de owes hg ace 
end The time for the onset of crazing for different geome tte # was general- yur om or heat-r : st 
stresses 1s found by obsers Ing the specimen periodi- crade cen standard Hmm 8 and r wae tape 
cally and noting the extent of the crazing. Some of tensile specimens and an accelerated weather 
the details of the test procedure are as follows. The | c..cimen were taken from each disk. An | 
specimens being tested were observed with the un- sania! of control specimens were cut out of 
aided eye, and, in addition, in the tests on the sheet from a location adjacent to the piece used 
representative samples nearly all specimens were also r Nanaia Sat iad 
observed with a 20-power Brinell microscope (fig. 5 While it was desired to obtain the same amo 
Observations of the extent and nature of crazing of stretch to within about 5 percent on the tl 
were made several times on the first day the load was formed pieces of a sample, the actual variatior 
> Sn oes a through the — day, | clongation between disks was greater than 
amaaied ee en eee amount, except for sample Pla. The 
‘ > € ~, © CTaZ * . c¢ - ; ° 
mum cross section was measured with a paper scale individual disks, based on the formula, eq | 
to the nearest millimeter; this corresponds to an 
accuracy of better than one percent when converted 
to stress. The observed data on the extent of the 
crazing were converted into stress values and the NBS sampk meetes NBS sampk 
latter were plotted against log time 


follows: 


3.5. Accelerated Weathering Test Li 


L2d 


The accelerated weathering test employed was the 
on . In these references the reduction in thickness was reported, not t! 
sunlamp-fog chamber type. This test was made in | the elongations were calculated by using eq (1 
“ It was interesting to note that the elongation at failure was 


accordance with Method 6021, Federal Specification | oracr orabout 20 percent by the slight stretching 
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Next, it was observed that the formed specimens 
figure 6 (C, D, E, F), had a laminar 
material 


4. Results and Discussion 


as shown in 
fracture. The laminar fracture of formed 
indicates that the segments of the poly mer molecules 


4.1. Standard Tensile Tests 


preferred orientation in. the plane of the 
propagation on planes 


standard tensile tests on the have a 
sheet, thus favoring fracture 
plane of the sheet In the 
of the polymet 


esults of the 
nd the unformed portions of the four repre- 
nearly parallel to the 
unformed material the segments 
to be randomly oriented 


readily seen from figure 6 


samples are shown in table 1. Figure 6 
s the appearance ol the fractures on broken 
ns of formed and unformed material molecule are assumed 
nsile strength and secant modulus of elas- Although it is not 
the four samples of polymethyl methacry- | on the fracture surface of each specimen there was 
e unaffected by stretch-forming to about | a perpendicularly 


nt elongation. The elongation at failure, 
increased by 


small mirror-like area oriented 


r xploratory tests, Was gre atly 
from approximately 10 percent to about 
nt 
train at the threshold of crazing also was 
greatly as a result of the forming. In fact, 
ples L2d, Pla, and P2a, at least half of the 
ns showed no stress-crazing up to rupture 
on other specimens was very light or was 


| only at accidental finger marks 


b. Discussion of Fracture Behavior 
of interest to consider the fracture behavior 
fracture mechanism in the formed and un- 
material. First it was noted that while the 
ns of the unformed material commonly failed 


0 10-percent elongation with the fracture | Ficure 6. Effect of biaria 


elongation on the 
mately flat and normal to the tensile load line, ee 
A in figure 6, occasionally a specimen | ,,.4. 7° etention bt werent 


imple I 


elongation, sometimes | elongation 77 percer 
wercent sarnple xi. forme 49 perce 


stretch-forming to about 50 percent 


¢ &pe mens 


n by 
d a much greater 
inied by an oblique fracture (B in fig. 6 
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to the tensile load line. On specimens of unformed 
material, the mirror area was diffusely bounded; 
on the formed, the boundary was sharp. The 
mirror-like surface was of the order of a few square 
millimeters in area and was found extending inward 
from the boundary of a cross section 

On unformed material this area was located at 
the corner of the cross section for most of the speci- 
mens in the other instances if was found either 
extending inward from the cast edge or from the 


machined edge At corner locations the murror 
was roughly a quadrant of a cirele, and at the edges 
it was semicircular 

On specimens ol formed material the mirror area 
was low ated at a corner on almost all of the speci- 
mens Howe ver, on the other specimens the mirrors 
were located only at the machined edge At corner 
locations the mirror area was almost a right triangle 
with the hypotenuse a convex curve instead of a 
straight line: the area extended much farther along 
the east edge than along the machined edge When 
located on a machined edge, the mirror area appeared 
to be somewhat less than a semiellipse, with the 
major axis normal to the edge In testing a formed 
specimen, usually some edge cracks appeared normal 
to the edge after the material had been strained con- 
siderably These cracks became larger ais the speci- 
men was stretched further, and at failure the frac- 
ture appeared to go through one of them. In figure 
6 two such cracks are evident on the left edge of 
specimen C just above the identifying letter. Such 
cracks reflect light similarly to large crazing cracks, 
indicating a mirror-like surface. Also, the large 
edge cracks extend farther in the direction of the 
cast surface than in the thickness direction, as in the 
cases of the mirror areas on the fracture surface 
For this reason it is quite plausible to expect the 
fracture in specimens of the formed material to 
start at such edge cracks 

From the fracture behavior of specimens of un- 
formed material, discussed below, it is logical to 
suppose that the fracture in such material also 
begins in the mirror area and that in this material 
this area is an extension of a crazing crack. 

The experimental evidence, that suggested that 
the fracture started at the mirror area is as follows 
The fractures were examined on a large number of 
unformed tensile specimens which were solvent- 
crazed and then broken. The tensile tests were done 
in connection with another phase of this investiga- 
tion on crazing [5]. The specimens, of the standard 
tensile type. were wetted with benzene while under 
a load; the solvent was applied to the central 4- by 
2-in. portion of one face of the reduced part of the 
specimen. On all the specimens inspected it was 
noted that a semicircular mirror-like area was present 
on the fracture surface; this area was located with 
its center at or near the solvent-crazed surface. It 
is plausible to suppose that this mirror area is an 


’ In this report these edge cracks have mev t arbitrarily been considered as 
listinct from crazing cracks; also any fine « “” on machined edges was dis 
regarded The reason is that suct racks and crazing are dependent on the 
machining of the specimens, and could pos y be minimized or caused to begin 
at higher strains by varying the machining technique or by properly annealing 
the specimens 


extension of a solvent-craze crack, and hy 
fracture is initiated at such a crack. This | 
sis was strengthened by comparing the locati 
fractures on a number of these specimens wit} 
graphs of the solvent-crazed specimens tal 
to breaking In all cases, the fracture was 
pass through a crazing crack. Visual evi: 
justify further this supposition was obtair 
different portion of the investigation, to be d 
in detail in a separate report. In the latt 
experiments were made on several different ¢; 
methyl methacrylate sheets of viscosity 
molecular weights ranging from 90,000 to 3 
In stress-solvent crazing tests on low-m« 
weight material, one or a few large crazing 
developed, and the specimens were seen to 
the rapid growth of one craze crack. It sh 
noted that solvent crazing reduces the tensils 
of the specimens [5] 

Because the evidence on the solvent-crazed 
mens strongly suggests fracture propagating 
ward from the mirror area that originates at 
ing crack, and because similar mirror areas extet 
in from the edge are found in standard tensile sp 
mens that are not solvent-crazed, it seems reasona 
to conclude, as previously noted, that failure beg 
at the surface of the specimen and that the mecha 
isms of fracture and crazing are closely relat 
Indeed, one might go further and sav that in 
specimens of unformed material that are not solver 
crazed the fracture first starts at a stress-crazing 
crack. In the formed specimens, which frequent 
do not exhibit stress crazing, the fracture of 
specimens is delayed, that is, occurs at a much hig! 
strain than in the unformed material. Furtherm 
the true stress is probably higher at or near fa 
in the formed material than in the unformed as 
actually found in a few tests.* 

The difference in appearance between the n 
part of the fracture surface and the rougher port 
mav be associated with a low velocity of fract 
propagation at the former and a higher velo: 
the rougher portion. Such an explanation 
fracture behavior of glass is discussed by Mor 
in his monograph on glass 

The fact that on the broken specimens of | 
material the mirror area had a smaller dimens 
the direction perpendicular to than parallel 
plane of the sheet, or laminas, suggests tl 
rate of crack growth, and perhaps the subse« 
high-speed fracture, too, is slower across the | 
than parallel to them 

Hsiao and Sauer [7], who studied cra 
specimens of polystyrene, present a different 
of the relation between crazing and fracture 
conclude that “the fracture cracks of the n 
are not the same as crazing cracks and that 
* In the exploratory work with specimens of a disk biaxially stretct 
cent elongation, the load-elongation graph was taken out to 50-percent 
3s low-magnification recording strain gage After the maximum 
tained, at about 4-percent strain, the load declined with increasing stra 
four-fifths of its maximum value and then remained almost consta 
percent strain, the limit of elongation of the strain gage; at this point t 
increasing slowly A3 the volume of the material remains practically 
the plastic range, the reduced cross-sectional area can be calculated fr 


and the true stress then derived. It was found that the true stress 
strain was about 10 percent greater than at the maximum boad 
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fracture is usually some flaw in the materia! 
ne of the crazing openings.”’ However, it 
reasonable to expect that, the 
crack produces a stress concentration at 
the subsequent fracture would be initiated 
\ microscopic examina- 
by 


mn since 


ex of the crack 
fracture surfaces 
8] on 


such as carried out 


coworkers several materials, 1s 


this study may clarify the situation 


4.2. Stress-Solvent Crazing Tests 
reshold data for the stress-solvent 
ests on the 50-percent biaxially stretched 
shown in table 2. Typical specimens are 
figure 7. The threshold was determined 
using two criteria For the first, called 
\, the threshold was taken as the maximum 
ow which there was no regular distribution 
eracks visible the unaided eye, the 
cracks being disregarded This was the 
terion used in the exploratory tests For the 
B, the threshold stress was taken 
crazing 


stress 


to 


riterion 
maximum which 
re vy isible to the unaided eve 
ght be expected, the average threshold 
stress obtained by criterion B is slightly less 
criterion A. The principal results as ob- 
v the two methods are in agreement, how- 
d are as follows 
verage threshold crazing stress for general- 
grade polymethyl methacrylate, crazed 
nzene by the procedure described previously, 
2.000 |b/in.? for unformed, and about 3,400 
) lb/in.? for the 50-pers ent biaxially stretched 
The values for the heat- 
gcrade 3.000 Ib/in and 5,000 
lb/in This represents an 
ment 80 percent for each 


below no 


stress 


corresponding 
are about 
respectively. 


from 70 


) 


to 


of 


or stress-solvent crazed s 


T) 


1< 


ipplied toa No 
i marked glass 


Was ult 
om pper 


al b44- by 3-in } f 
er load, was stroked wi 


and then were rem< 


i bri 


th the The specimens were under 
| testing machine; after 4 to 
examined for crazing thresh The crazing was done in a con 
nperature of 23° C and 530 percent 
for 1 week in this atmosphere 
{ the sample 


ved fron f 
were ‘lds 
phere room operating at a te relia 

iit 


ifter conditioning the s 


pet 
formed disks, 1 from « 


the tapered tensile specimen, 


Not only did the formed specimens exhibit higher 
threshold stresses than the unformed, but also there 
was a tendency for the crazing cracks to be somewhat 
finer and more closely spaced on the formed speci- 
mens 


4.3. Long-Time Tensile Test 
a. Threshold stress-crazing data 


crazing at 1, 
plots of threshold 


stress 


Values of threshold stress for 
10, and 100 hr, derived from 
stress versus log time, are given in table 

The plots of threshold-crazing Stress versus log 
approximately linear For 


time were in general 
most specimens the extent of crazing was recorded 
separately for each face of the specimen, because the 
crazing often progressed more rapidly on one face 
An examination of the plotted 
data showed no consistent behavior of the masked 
relative to the unmasked face Any consistent 
difference in the extent of crazing on the two faces 
was accordingly assumed to be caused by a slight 
misalinement of the specimen. In such a 
single straight line was fitted to the data for the two 
} were taken from this 


than on the other 


Cases 


and the values in table : 


faces, 

line 
The slopes of the threshold stress versus log time 

plots appeared to be about the same for all materials 


and test conditions. The effect of such factors as 
forming, relative humidity, sample, and grade on 
the threshold stress is indicated best by examining 
the 100-hr unaided-eve values, as these are the most 
numerous and These values, taken from 
table together with values of the ratio of the 
threshold stresses, S., for the formed and unformed 
specimens of each sample are listed in table 4. These 
data indicate that forming i biaxial strain of 
about 50 percent increases the threshold-crazing 


precise 


, 
>, 


1o 


per 


or a 
rhe 

Wf sheet-t 

wh 


ssib 
That poir 
sible with the unaided ¢ 
1 B: That point below whi 


4 Criterior low 
cracks Vv 
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* One 
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t craze under t 
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f all samples about 40 to 50 percent From loading specimens The specimens were photo- 
is also evident that the threshold stress is | graphed while under load as the finest crazing cracks 
} to 50 percent higher for the heat-resistant | usually were not visible on removal of the load as 

the general-purpose-grade samples. The | noted by other workers 
not precise enough to determine with cer- The threshold-crazing-stress values give an incom- 
nv difference in threshold-crazing stress due plete picture of the effect of forming and other 
ve humidity or to material of a given grade. | variables on the crazing behavior of the materials. 
Thus, the crazing on the formed specimens, where it 
b. Appearance of Specimens occurs, Is usually finer than that on the unformed 
specimens. Also, although the threshold stresses at 
95- and 50-percent relative humidities did not differ 
appreciably in general, the nature of the crazing at 
Threshold-crazing stresses of polymethyl methacryl these two humidities was markedly different The 
Poet petite, A rg var cracks at 95-percent relative humidity appeared finer 
and more closely spaced, and were almost always 
noticeably shorter than those at the lower relative 

humidity 

For most specimens the lengths of the longest 
cracks in the vicinity of the minimum section were 
measured with the Brinell microscope. A few cracks 
of exceptional length, appearing to be two or more 
cracks joined together or to be initiated at very fine, 
long scratches, were disregarded The cracks meas- 
ured on the unformed specimens after 100 hr at 
50-percent relative humidity were from 0.4 to 0.7 mm 
long; at 95-percent relative humidity they were in 


? 


general about 0.1 to 0.3 mm long The formed 
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8, taken near the end of the testing period, 
es the appearance of the long-time tensile- 


sampies ajle LU) 





FIGURE Ss pec 
am ple P2a, 
informed, afte 
fensile loading 
and 5 000 } , 


the minimum sect 


formed and wunfi 


mens, respective 


specimens, in the few instances where the data were | formed and unformed specimens of samples | 
available (all on heat-resistant material), had corre- | Pla at the high humidity was about the sar 
sponding crack lengths of close to 0.2 mm at the | stress at the minimum cross section was the same { 
50 percent and 0.1 mm at the 95-percent relative | all specimens of a sample. Sample Pla is t] 
humidity one of the four materials on which it was not poss 
It should be noted that the formed and unformed | to obtain values of threshold stress for stress craz 
specimens cannot be compared on the basis of the | at the high relative humidity; at the loads use« 
above crack-length data because the cracks on the material creeps before crazing can begin 
two se ts of spec imens had started at different times and 
the length of a crack appare ntly de pe nds on the “‘crack 
lifetime’, that is, the time elapsed under load after 
the crack appeared Nevertheless, when the crack 
lifetime is taken into account the cracks on the un- The effects of biaxial stretch-forming 
formed specimens seem to grow more rapidly than on | crazing behavior of polymethyl methacrylat 
the formed specimens even when the stress is | may be explained qualitatively on a molecu 
appreciably higher on the latter. For example, on | as follows. In the unformed state the 
one formed specimen ol sample P2a at 50-percent molecules are assumed coiled in an approx 
relative humidity, cracks were first observed at the | spherical shape; the chain segments hav: 
minimum cross section (6,000-lb/in.* stress) at about | ferred orientation. In the formed state tl 
50 hr; and after 50 hr of crack lifetime, the longest | cules should be somewhat uncoiled and in : 
cracks in this area were close to 0.15 mm in length. | disk-like shape with the chain segments 
On the corresponding unformed specimens with only | predominantly in the plane of the materia 
5,000-lb/in.? stress at the minimum, the cracks were | following mechanism of crazing, somewhat 
0.35 to 0.4 mm long for the same crack lifetime. The | to that proposed by Maxwell and Rahm [9 
difference in length is not as noticeable on specimens | lated. The crazing is assumed to start at th 
at 95-percent relative humidity at submicroscopic flaws or weak points. Suc 


4.4. Discussion of Mechanism of Crazing 


Another effect of the high humidity was to increase | points may be submicroscopic regions in wl 
the rate of cre« p markedly at the stresses used in the chance, the polymer chain segments are ¢ 
long-time tensile tests. Some of the specimens at | normal to the applied tensile stress. With s 
the 95-percent relative humidity “necked down” at | stress, a separation between portions of 
the minimum section near the end of the testing | chains occurs; a stress concentration exists 
period. The formed specimens of the heat-resistant | apex of the crack, and the latter grows 
materials necked down sooner than the corresponding | reaches a region in which the polymer cha 
unformed specimens. However, it should be noted | ments are oriented approximately in the d 
that these formed specimens had a 20-percent-higher of the tensile stress. The crack either does n 


stress than the unformed. The necking down of the ' or grows slowly unless the tensile stress is 
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Subsequent crack growth may involve 
‘f primary valence bonds, especially if the 
elatively high, of the order of the tensile 


cess Of biaxial stretch-forming, by orient- 
hain segments in the plane of the sheet, 
he proportion and size of the weak, normally 
regions and increases the regions of pre- 
ly parallel orientation. Stated differently, 
tch-forming may be introduce 
’ planes in the plane of the sheet. This 
yn or introduction of cleavage planes greatly 
the development and growth of crazing 
Thus, as neted previously in the long time 
aading tests, the crazing cracks, after be- 
visible, grew more slowly on formed as 
d to the corresponding unformed specimens 
ard to stress-solvent crazing, mechanisms 
en suggested by various authors [10, 11, 
h, while differing in some aspects, include 
tor the concept of the solvent acting as a 
er. By using this concept the mechanism 
be modified 


said to 


d above for stress crazing may 
le the influence of solvents as follows The 
molecules penetrating the surface of the 
tend to surround portions of the polymer 
ind reduce the forces required to separate 
Because of this weakening influence of the 
molecules at a surface flaw. such as a reqion 
al orientation of the polymer chains, the 
meentration that can be withstood is reduced 
ny crack develops at a lower applied stress 
the absence of solvent The solvent mole- 
capillarity probably fill the crack as it grows 
tinue to exert a weakening influence at the 
In this connection, it has been suggested by 
ns, Baker, and Howard [12] that another 
ng influence at the apex of a crack is the film 
ng pressure of the crazing liquid. 
effect of forming on stress-solvent crazing 
De expected to be similar to that for stress 
The reduction in the number and size of 
vions of normal orientation and increase in the 
of parallel orientation should result in higher 


PaABLe 5 


threshold stresses for the formed material. Also, 
for formed material as for unformed, the crazing 
stress should be lower in the presence of than in the 
absence of solvent, owing to the weakening influence 


of the solvent 


4.5. Accelerated Weathering Tests 


The results of the 480-hr sunlamp-fog accelerated 
weathering tests are shown in table 5 

The light-transmission values for all the materials, 
both unformed and formed, are 92.0 + 0.1 percent 
initially and 92.3 +0.3 percent after weathering 
While the transmission values are slightly higher 
after weathering the materials, the individual differ- 
ences, which do not exceed 0.5 percent, are considered 
within the experimental error of measurements made 
at different times. The haze values are approxl- 
mately 0.5 +0.2 percent for all materials both 
formed and unformed and before and after weather- 
ing. In this connection, the specimens were in- 
spected visually after the weathering test, but no 
crazing Was observed on any of them 

For all samples, the shrinkage of the unformed or 
control specimens was very slight, averaging 0.05 to 
0.1 percent, which values are of the order of magni- 
tude of the standard error of the shrinkage values 
The formed specimens of the heat-resistant grade 
samples shrank only 0.2 percent; however, similar 
specimens of the general-purpose grade samples 
shrank somewhat more, the values being 1 and 2 
percent for Plexiglas I-A and Lucite HC201, respec- 
tively 

Most of the specimens were slightly warped afte 
the weathering test. In general, the formed speci- 
mens were more warped than the unformed, particu- 
larlv for the general-purpose-grade materials 

The greater dimensional changes for formed pieces 
of the general-purpose-grade samples as compared 
to the heat-resistant grade are not surprising for the 
following reason. The specimen temperature in the 
test 7 which is not far below the 
second-order transition temperature of the general- 
namely, 75° to 80° C; the 


is about 60 5° (' 


purpose-grade samples, 





temperature for the heat-resistant 


to U5 CC 


( orresponding 
crade is about 94 


4.6. Possible Applications of Stretch-Forming 


crazing properties of 
by moderate (50 


The 


acrvlic 


improvement in the 
plastic produ ed 
biaxial stretch-forming suggests that formed 


sheet 
percent 
enclosures made from prestretched flat sheets may 


have superior crazing and strength properties com- 
pared to enclosures formed from normal sheets. In 


an astrodome as prepared normally, as, for example, 
by vacuum drawing, there is a maximum amount of 
stretch and craze resistance at the apex and negligible 
stretch and minimum resistance at the rim 
The use of prestret hed sheet would hence improve 
the craze the rim where 
contact with crazing liquids is quite likely 

As an alternative to using prestretched material to 
achieve improved craze resistance at the edge of 
acrylic enclosures, there is the possibility of preparing 
an enclosure larger and more deeply drawn than 
required, and then using only the central portion of 
the formed piece 

In view of the considerable « hanges in the phy sical 
properties of polymethyl methacrylate 
resulting from biaxial stretching to about 50 percent, 
it is of interest for both practical and fundamental 
reasons to learn the effects of greater stretching on 
this material. Experiments with this objective are 
in progress on material stretched up to about 150 


craze 


resistance, especially at 


cast sheet 


percent 

The effects of biaxial stretching on other physical 
properties, such as creep, abrasion resistance, natural 
weathering while under impact resistance, 
fracture under bullet impact, and crazing with 
various solvents, should also be determined 


load 


5. Summary 


Some of the effects on polymethyl methacrylate of 
biaxial stretch-forming to about 50 percent are as 
follows 

1. The strain at the onset of crazing in the standard 
tensile test is greatly increased; in fact, most speci- 
mens showed no crazing. The tensile strength and 
secant modulus of elasticity are unaffected. The 
elongation at failure is increased from approximately 
10 percent to about 60 percent 

2. The threshold stress for stress-solvent crazing 
with benzene is increased about 70 to 80 percent for 
both general-purpose and heat-resistant grades 

3. In long-time tensile tests of up to 7 days’ dura- 
tion (a) the threshold stress for stress crazing is 
increased about 40 to 50 percent for both grades of 
material at both 50-and 95-percent relative humidi- 
ties, and (b) the crazing cracks produced are some- 
what finer and appear to grow in length more slowly. 


These values were derived from volume-temperature measurements made 


on these samples over the range 2 to 110" ¢ s mercury dilatometer was used 





The effect of biaxial stretch-forming 
physical properties, as well as the effects 
degrees of stretching, should be investigat: 
practical and fundamental reasons 

The considerable increase in the 
failure and in the stress crazing and str 
crazing threshold of polymethyl methac: 
of moderate 50 percent) 
suggests that formed enclosures made 
stretched flat sheets mav have greatly 
crazing resistance and strength properties 
to the extent that acrvlic glazing need no 
inated 


elor 


result strete 
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Mass Spectra of Fluorocarbons’ 
Fred L. Mohler, Vernon H. Dibeler, and R. M. Reese 


orocarbons 
are present 
I methane 
neXxane 
lifferent 


four olefins 


l. Introduction 
wing interest in fluorocarbon chemistry 


ology makes it worth while to obtain mass 
perfluorohydrocarbons as a basis for chem- 
sis of unknown mixtures of fluorocarbons 
theoretical int how 
ffer from the 


SS Spectrometry sectlion of 


s of erest to see these 


hydrocarbon mass spectra 
the National 
if Standards has for several years been col- 
ass spectra of fluorocarbons as compounds 
and the paper 
on of spectra obtained to date A previous 
includes data on five perfluorohydro- 
The present includes 
ids Other laboratories have obtained mass 
i a variety of fluorinated hydrocarbons, but 
compilation of spectra has 


ivailable, present vives a 


summary these 


matic been 
lete spectra of all the fluorocarbons included 
iper have been published in the API Catalog 
Spectral Data [2]. When the mass spec- 
a fluorocarbon is tabulated in the conven- 
inner in the order of increasing molecular 
f the ions, ions with different numbers of 
atoms and of fluorine atoms become mixed 
For that reason, spectra 
ented here in order of the number of carbon 
nd the number of fluorine atoms in each ion 
convenient for the purpose of interpreting 
tra in terms of molecular structure and for 
son with hydrocarbon spectra. For pur- 
chemical analysis the conventional tables 
ourse, preferable 


nfusing disorder 


2. Experimenta! Details 


of the mass spectra were obtained with a type 
solidated mass spectrometer equipped with a 
anometer [3] to measure the pressure in the 

was presented tu itt irgh onference or r tica 


ctrochemi 
ts indi 


inlet reservoil \n voltage of 70 volts was 


ised, and the metastable ions were suppressed Kor 


lOnIZINEG 


the heavier compounds thre suppressor voltage “Was 


adjusted to give a resolving of about 350 
The temperature ol the chambet 
250°C. The heaviest compounds require a magneti 
field in excess of the rated capacity ol the magnet 
but there was no evidence of overheating in the time 
of the 


Miwass 


powell 


ionization was 


some 
type 102 
instrument Was In- 


required to scan the spectrum mass 
with a 


103 


were obtained 
before the 
stalled. There are only minor differences in spectra 
obtained with the two types of instrument 

Weare indebted to Dr Creorge H Cady, | niversity 
of Washington: The M. W. Kelloge Co.: Minnesota 
Mining & Manufacturing Co.: the Jackson Labora 
tories of the E. 1. du Pont de Nemours Co., Ine.: and 
to others acknowledged later for making available 
the various fluorocarbons. Most of these compounds 
are available only in small experimental batches 
The purity was in all cases adequate for a reliable 
description of the spec tra Trac e amounts of mole 
cules containing one H or one Cl atom were som 


spectra spe 


trometer type 


times found, but, except where noted, there was no 
evidence of heavier fluorocarbons as impurities. The 
API tables include mass peaks as small as 0.01 pet 
cent of the maximum peak. For brevity, we omit 
many of the small peaks in the tables to give a con- 
cise presentation of the major features 

There are no difficulties in running fluorocarbons in 
the mass spectrometer, although special precautions 
were taken in handling the iso-C,Fs because of the 
reported toxicity of this compound [7]. The vapor 
pressures are high and are of the order of magnitude 


| of the vapor pressures of the hydrogen analogs 


Thev are inert and are not strongly adsorbed, so they 
are quickly pumped out, with no serious memory 
effects in spite of the high molecular weight 


3. Results 


Table 1 lists the larger ion peaks in the mass spec- 
tra of normal perfluoroparaffins from methane to n- 





TABLE 1 


heptane. Ions are listed in the order of the number 
of carbon atoms and the number of fluorine atoms in 
each ion, and intensities are relative to CF? taken as 
100. The last row gives the sensitivity (current per 
unit pressure) for the CFy ion relative to the sensi- 
tivity of normal butane at mass 43 

In each of these spectra by far the largest peak is 
CF;, and the molecule ion peak is very small or 
zero. 
atoms the maximum peak is C,,F,,. 
molecules heavier than propane, 


In each group of peaks containing n carbon 
, except that for 
when one carbon 
atom is removed from the molecule the maximum 


peak is C,F,,_;. Thus in C,F,, the largest peak con- 
taining five carbon atoms is C;Fe. In C;F,, the 
largest peak containing six carbon atoms is C,F¥%, 
but the largest peak with five carbon atoms is C,Fy 
The size of these maximum peaks decreases progres- 
sively with increasing molecular weight, and the 
spectra of the heavier compounds are very similar 
except near the heavy end, and these distinctive 
peaks are small. In normal hydrocarbon paraffins 
the maximum peaks are also C,,H,,4,, except that the 
molecule ion peak is larger. The largest of these 
peaks is C,H, or C,H. 

In these fluorocarbons the sensitivity for the CF; 
ion progressively increases with increasing molecular 
weight. ‘Lhe total ionization (the sum of the mass 
peaks times the sensitivity much the 
same rate. 

Doubly charged ions are observed as follows: In 
CF, the intensities of CF}* and of CF;* are 3.9 
and 1.25, respectively, relative to CF; in C,F, the 
intensity of C,F7* is 0.28, and in C;F, the ion 
C,;F3* is 0.09. Thus, molecules of formula C,F, 
give ions C,F**, and in molecules where n is an 
even number a singly charged ion will mask the 
doubly charged ion. In spectra obtained with the 


increases at 


n+2 


Vass spect aot norma 


pe flluoroparaffins 


type 102 mass spectrometer some metastabl 
tion peaks were observed, but they wer 
small [1]. 

There are only two examples of branched flu 
carbons, perfluoroisopentane and isohexane. Ta 
2 shows the spectra and, for comparison, the spect 
of the normal compounds from table 1. The mos 
conspicuous difference in both the pentanes and thy 
hexanes is found in the relative intensities of 
C;,F; ion. It appears as if the isopropyl rad 
CF(CF;). did not contribute at all to this peak, al 
that the normal propyl ion CF;-CF,-CF; account 
for the observed ion current. In isopentane, CF 
is nearly zero, and in n-pentane and _isohexai 
each of which contain one n-propyl radical, 1) 
relative intensity is about 9, whereas in n-hexan 
which consists of two n propyl radicals, the relatiy 
intensity is twice as great. C.F? is also 
twice as great in the normal compounds containing 
two terminal ethyl radicals as in the iso-compou 
containing one terminal ethyl! radical. 

Table 3 gives mass spectra of three pure cy 
and methyleyclohexane. The spectra are \ 
distinctive with few common characteristics. T! 
molecule ions are small, and, with two minor exe: 
tions, the peaks C,F,,_, are the largest peaks for 
greater than 2. In the pure cyclics, production 
CF; involves a rearrangement of atoms, but it is 3 
large peak. The ion C,F; also involves a rearrange- 
ment of atoms in these compounds. No dou! 
charged ions are observed in these spectra 

Table 4 lists mass spectra of four perfluoroole! 
The spectra are distinctive, with no obvious comn 
characteristics, except that the molecule ions ar 
fairly large and CF* is large. CF is small 1 
C,F,, where a rearrangement of atoms is involved 
the ionization process. The two butenes have ver 


rougi 


344 








different from the 
The relative 
‘+ and CF; 
spectra C.F 
bon peak } Ik of 


clistinetive spectra and both 
evclo-C, F riven 
of CF C;k 
different in the thre 
a relatively large doubly charged 
intensity 2.17, and C IK 
l-butene C,F ol 
Table 5 lists mass spectra of three C.F 
1.3-Butadien CF Ck Ck CF cvclobutene 
CF,—CF=CF—CF,; and 2-butyne, CF,—C =C—CF 
In thes 
The relative 
distinctively 
Table 6 


a trievelu 


isome! n tat ; 


intensities are 


rives 


strikingly 


has intensitv 0.16. In 


iIntensityv U is observed 


compounds C,F? is the maximum peak 
‘*, and CF 
ee compounds 
spectra of a dicevelic and 
The dievelic compound is 
completely fluorinated 1-methylnaphthalene, C,,F 

prepared by E. T. McBee of Purdue University 
This molecular weight 512 is the 
heaviest compound ever run on the Bureau’s mass 


ascribed to & F 


heaviel 


intensities of ‘ Ck are 
different in thi 
Mass 


irk ules 
fluorocarbon 


compourn of 


wo small 
that ther ~ 


spec.romete! DCaks 


and C,,F 
Impurity 


indicate some 


The maximum peak of each 


present I 
formula C,F for n greate! 


group has the or 
The ions C.F? and C.F? involve rearrange- 
im the There are 


included in table 6 that come either 


Cul bon 
than 4 
ment small 


Irom 


ionization also many 
pe aks not 
rearrangement or [rom impurities 

The C.F). compound was furnished by W. T 
Miller ot Cornell | niversity and Was described by 
Prober and Miller in 1949 [4 When hexafluoro- 
butadiene is heated, one of the products is a saturated 
fluorocarbon of formula C.F It has a freezing 
point of 40° C boiling point of 80° C. They 
that ‘“‘preliminary X-ray diffraction studies are 
consistent with the ; 


and a 
State 


following) structure 


The mass spectrum of table 6 can also be described 
this structure While the 
CF? and C.F? involve rearrangements in the 
ionization process, if the structure, other 
spectra show that rearrangements to give these ions 
are very common in fluorocarbon spectra. It is 
perhaps significant that C,F? is the largest C, ion, 
C,Fy the largest C, ion, and C.F? is fairly large 
The spectrum shows a very uniform distribution of 
intensity from CF* to C.F The largest peaks in 
each carbon group are C,F),,_; beyond n=4 


as consistent with ions 


this is 


4. Conclusions 


The fluorocarbons included in this investigation 
cover a wide variety of compounds, but, with the 
exception of the normal perfluoroparaffins, there are 
not enough compounds of any chemical class to 
permit broad generalizations. The normal perfluoro- 
paraffins show a marked regularity in their spectra 
that is quite different from the regularities found in 


the hydrocarbon analogs CF is by far 
abundant ion and accounts for half or m 
total ionization Anothe1 
the appearance potential of this ion is 
14.3 to 14.4 ev) in C.F,, C3F;, C.F 

The value riven in the reference tor ( 
As produc tion of CFF is the pre 
is significant that tl 
Increase 


research ipl s 


correct 
mnization 
ity, or current, pet 
ively with increasing molecular weight 

In the ionization valence 
removed, and this weakens the ionized bor 
either the mole 
was initially 
themselves bet 


proce ss, if 
unit pressure 


process a 


are two possibilities, l 
ciates at the point where it 
2) the electrons rearrange 
is time for dissociation and the disso« lat 
at some other bond than at the point initial 
On the first hypothesis, the probability of 
CF? should be the 
heavier than ethane, for 
minal C—C 
seen that the ionization potential remains t 


and the progressive 


same for all normal co 
Ih eae h Cause one of 
bonds has to be ionized 
This is not the case, 
the probability of ionization indicates that 
the molecule is initially ionized there is a | 
ability of breaking a terminal C—C bor 
CF Large molecules with many vale! 
afford a larger target area than small 
Wallenstein, Wahrhaftig, and Eyring 
vanced this second hypothesis concerning 
tion process to explain the relative intens 
hvdrocarbon mass spectra, but in hydrox 
radicals other than CH, are preferentially 

It has been noted that in normal perfluoro) 
ions of formula - Kr rive the largest peat 
carbon group, except that in fragment ions 
removed from the molecule C,F 
peak. In hydrocarbons, C,H, 
largest, with rare exceptions. In the ey 
ions C I, tend to be the largest for 
than 2. In the dicyelic, C,,F., the ions C,] 
largest, and in the tricyelic, C.F,, 
C,F; tend to be largest. In this respect 
of fluorocarbons are simpler than hvdroe 
In hydrocarbons there is a close resemblance 
patterns of the 1l-olefin and the 
There is no such resemblance between ©) 
and perfluoro-1-butene 


carbon 


largest 


ions ol 


isomer! 


Identification and analysis of mixtures of px 
paraffins will be difficult because of th 
complete absence of molecule and 
intensity of all the heavier ions in these 
Unsaturated molecules and cyclics have vi 
tinctive spectra and possibilities of analysis 
favorable when spectra of more compounds 
available. 


1oOns 


As there are gaps of severai mass units 
fluorocarbon peaks, impurities other than 
carbons are usually conspicuous. A commo! 
impurity will be molecules that are not co1 
fludrinated and contain one or two H atoms 
previous paper on mass spectra of fluoroca: 
includes a spectrum of cyclo-C FH, and it 
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peaks contamming an H atom correspond 


relative intensity to large peaks in the 
spectrum, with H substituted for one of 
is We have qualitative evidence based 
n impurities rather than pure compounds 
s a quite general relation and also holds 
atoms are substituted for 2 F atoms in a 
Thus, if a trace of H is present, peaks 
H will not be adjacent to the large fluoro- 
iks but will be 18 mass units less than the 
S With traces of chlorine in the fluoro- 
larger peaks containing Cl will be 16 and 
nits heavier than the larger fluorocarbon 


nv HF or F, is formed in a reaction this 
glass to give SiF,. The Sif? ion of mass 
ar the largest peak ip the Sif, spectrum 


l 


ll be a common impurity peak 
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A Method of Computing Exact Inverses of Matrices 
With Integer Coefficients’ 


J. Barkley Rosser’ 


he problem of computing the exact inverse of a matrix A with integer co- 


} ry ith 
efficients is completely solved by solving exactly the simultaneous equations A 


¥, in whic! 
both z and y are variable vectors This solution can be carried ¢« 


ut bY any one of numerous 
well-known procedures, resulting in expressions for the components of x as linear combina 
tions of the components of y. The coefficients of these linear combinations are just the 
components of A~! because we have 1-ly 

In actual practice, if the order of A is at all large, the exact components of A~! will be 
fractions whose numerator and denominator each have a large number of digits, and the 
isual methods of solution become extremely laborious due to the necessity for carrying an 
even larger number of significant digits throughout most of the computatior In the method 
presented herein, the number of significant digits involved builds up gradually, and only the 
final stages of the computation involve a large number of digits. Moreover, the method 
can be readily adapted to use on IBM equipment, and so all but the 
many significant digits must be carried) can be readily mechanized 


final stages (in which 


1. Illustration of a Solution by Previous Methods 


Suppose we require the exact inverse of 


128 ISS 


To invert this, we write 


OLr, 
Li3r, 
1382, 


120r; + 102r, 


Mir 772, — 154234 1177,— 1312 1122 


Eliminating z, from the last five of these by use of eq 1 gives 
344z; 48637; + 2238 lr, — 350437, + LIS Il r6= 103y 
18577,+ 222r,— 1400r,+ 22847;— 36257%,= y 
9088r, + 343557, + 340072, — 329292, + 229692, = 157y 
166162,.+ 122317 63412, 52092, 53992, 


361527, + 583617 19312, 70197. + 118672. 191 y 


Che preparation of this paper was sponsored (in part) by the Office of Naval Research 
National Bureau of Standards and U niversity of California at Los Angeles 
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The multipliers that were used in eliminating r, are apparent from the coefficients of 
y’s on the right sides of the reduced equations 


We now use eq 2 to eliminate 2 from the four equations that follow it. The multip 
| i | 


we use for eq 2 are, respet tively, 4857, 1136, 2077, and 4519 We get 


OS2 22017 1694 18155 561 19027 5 03367 7 38264 6536, 


383 U29001 124 29629 1 72672y H530y, 
2674 730 08298 242 992900, 3 15704y.— 3268y 


0 56706 1586 61134 538 84100 : 6 SOSRSSy 6536y 


We now use eq o LO ¢ liminate z, from the three equations that follow it The multiplic 


we use for eq 3 are, respec tivelv, 3 68507, 1 67974, and 10 24546 We get 


LOLS 85427 75664 185 84797 36168 527 848906 85006 t SS720 40000, 


5 809046 41224y 24085 61L752y SO0OUS8S8 GS712y 


61 00170 O1012 $1112 09360y 
63889 07320y 10978 78064 y 13498 11452 


90 GO003 96128 160 74202 20652 


1434 17256 47920 2302 30153 46048 927 22912 23488 7 12009 O3008y 


9 47491 69752 66964 326564 SOUSS 68712, 


dal/ 


We now use eq 4 to eliminate z, from the two equations that follow it. The multiplic 


we use for eq 4 are, respe tively, 63 21634 and 997 39385 We get 


26343 37340 69713 22772 854 75613 37432 235402 17 65250 77603 O4880, 85 44010 53886 36424 
93 01742 24170 60256y 51 19808 38113 15408y, +95 64239 81708 51S884y;, 


171 10024 03515 99936y 


1 49345 83813 20139 OS89036 | 30525 23686 92201 70136 
71168 36351 59272, 807 77618 45306 22120y,— 573 85438 90251 


919 23447 01231 65744y 714 


We now use eq 5 to eliminate z; from the equation that follows it. The multiplier we uss 
for eq 5 18 13 52291 10622 We get 
1 78080 76167 83302 24294 95529 751042 646 84381 45668 50109 91756 S8544y 
1037 27853 59149 81148 61699 71LOOSy 146 95683 03519 19158 33166 05736, 
1234 46556 51847 48058 62839 54504y,,4+ 1293 36364 42399 49987 74131 11848, 


1368 83204 20859 69018 78357 06376y, 


Dividing this through by the common factor 191 28479 63417 03768, we get 


930 97185 49728 3 38157 46308 J; 5 42269 20056 IP 2 33660 40527y 6 45354 77403 44 6 76145 55311 


7 15598 97507, 


Substituting this back into eq 5, 4, 3, 2, and 1, we get 


97185 49728 | 88884 62180, 99087 63656y.—2 41312 82729y 


60486 42243y,+5 90937 51511ly,—2 67697 43115 


97185 497287, 5 46840 79248y 1 05392 26192y t 94982 OS180y 
3 20138 43300y,+— 5 11496 67628y 1LO65R S5132y 
97185 49728: 29934 7972 4 29235 11632y 13089 17764y 
3 68561 73708y,—2 89478 23484y | 30251 85404y 


97185 497282 1 95059 83764y, +5 01313 85144y,.+2 21002 58797, 
3 77781 95727y,—3 87417 30707y; + 4 19819 57079 


O7185 497282 | 68307 83936, $1475 04752y.+2 21089 12392y 


1 71153 69576y,+ 4 11243 40200y | 67986 90584 


From these results one can easily write down A 
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We hardly need to stress the fact that th procedure outlined above is not a really prac- 
le method to find .A™. In the present case, we did carry the computation through to the 
r end, just to show how unwieldy it becomes, but it required 38 hours of computing time 
trained professional computer. We might note further that the computation would 
been even more unwieldy and extensive if we had not deviated from a strict mechanical 
edure by removing common factors from our multipliers when combining equations 
iselvy, what we did do in this direction 1s described below 
In general, if we wish to eliminate z,; from two equations az br. , the 
hanical wav to proceed is to multiply the first equation by 6 and the second equation by 
and add However, in order to mitigate somewhat the frightful increase in the sizes of 
coefficients in the later stages of the process, we have in every such step of the present com- 
ation determined the greatest common factor. c. of @ and b, and have multiplied the first 
ition by 6/e and the second by —(a/c), and added. For example, in the case of eliminating 
om eq 5 and the equation following it, @ is a 20-digit number and 6 is a 21-digit number 
vever, the common factor, c, is a 10-digit number, so that our multipliers, 6/¢ and ale 
2- and 11-digit numbers, respectively. If we had used 6 and —a as multipliers, the coeffi- 
ts in eq 6 would each have had about nine more digits 
Similarly, in the multipliers used with eq 4, a six-digit factor was removed, and in the 
tipliers used with eq 3 a two-digit factor was removed. Without such removal of factors 
coefficients of eq 6 would each have had about 15 more digits, so that it is doubtless worth 
le to carry out such a determination of common factors Nevertheless, this portion ol 
computation can be quite a chore, especially when (as in the present case) the greatest 
nmon factor of a 20-digit number and a 21-digit number is required. Also, if determination 
common factors is carried out, mechanization of the process is more much difficult 
One could effect a further saving in the number of digits carried by extracting common 
tors from each of the equations derived in the course of the computation, instead of only 
mm eq 6. However, this entails a great increase in labor with only a moderate decrease in 
number of digits in the various equations. In the present case, if all possible common 
tors would be extracted from previous equations, one would still obtain as the equation 
responding to eq 6 an equation with a 23-digit coefficient for z,, and there would be much 


litional labor in the determination and extraction of common factors 
2. A New Method of Solution 


We treat the same matrix A as in the previous section. If G@ is the inverse of A, then GA 
the unit matrix. To find G, we find in succession B, C, D, &, and F with the properties 
hat BA has its first column the same as the unit matrix, CA has its first two columns the 
me as the unit matrix, DA has its first three columns the same as the unit matrix, and so on 
to GA, which has all its columns the same as the unit matrix 
Our method for finding B, C, , Gis a modification of the algorithm set forth in a 
revious note.’ In addition, we use the following well-known property of matrix multiplication 
Wand V are matrices with WA=V, then if we form W* and V* from W and V by performing 
same elementary transformation on the rows of each, we will have W*A=V™*. By an 
mentary transformation on the rows, we mean one of (a) Multiplying the ith row by a 
onstant a. (b) Interchanging the ith and jth rows c) Adding a@ times the ith row to the 
h row 
By a sequence of such transformations, we can reduce V to the unit matrix, and so the 
ime sequence of transformations performed on W will reduce it to A This fact is the basis 
various methods for computing A~'. The novelty in our method lies in the fact that we 
re able to use mainly transformation (c) with an occasional transformation (b) until the final 
age of the reduction, and also that we have a mechanical procedure for keeping the sizes of 


he numbers small until the final stages of the reduction 


Barkley Rosser, A note on the linear Diophantine equation, 





As indicated earlier, we carry out the reductions on V in a certain order. First we rex 
the first column of V to be the first column of the unit matrix. Then W has been red 
to the B mentioned earlier Then, restricting ourselves to transformations that leave 
first column of V unchanged, we reduce the second column to be the second column of the 
matrix. Then W has been reduced to the ( mentioned earlier. Proceeding in this way, col 
by column, we eventually reduce V to J and W to A 

To get started, let us take W to be the unit matrix, 7. Then Vis A. We first seek tr 
formations that will bring the first column of V to the desired form. So we temporarily igi 
all other columns of V, and consider only the first column, which is (1 This is to be redu 
to form (II 

I] 
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One can do this in many ways, but we follow the way that is proposed in the note referred 
in footnote 3, since this is quite mechanical but keeps the sizes of the numbers involved r 
sonably small. Specifically, we apply the elementary transformation (c) to those two rows 
containing the two numbers of maximum absolute value. Thus we first add 1 times tl! 
fourth row to the sixth, getting (III). Then we add —1 times the first row to the fourt 
getting (1\ Then we add —1 times the second row to the first, getting (V 


\ 


19 


34 


Clearly, if we were trying to reduce the sizes of the numbers as rapidly as possible, we would 
now add —2 times the first row to the second. However, this presupposes that good judg 
ment is to be applied at the various steps. One of the advantages of the procedure we a! 


describing is that it gives quite good results even when applied quite mechanically. TT: 


illustrate, we ignore the smart transformation, and proceed according to rule, adding 
times the ith row to the jth for the following values of i and j: 


then have 





now have a situation that occasionally arises, in which the largest number is at least twice 
In such case, our multiplier can conveniently be different from 

Thus we could conclude by adding 2 times the second row 
finally, by interchanging the first 


reat as the next largest 

1 used uniformly so far 
| times the second row to the fifth, and 
we could equally well continue mechanically, since exactly the 
to the first row, and then 


he first, then 
However 


second rows 
if we twice add 1 times the second row 


result would ensue 
interchange the first and second rows 


1 times the second row to the fifth, and, 


finally, 


(pplving the transformations listed to the unit matrix transforms it into B. which is 


| has the property that the first column of BA is the same as the first column of the unit 
trix. Moreover, the coefficients of B are quite small, which is why there is little increase 


the sizes of our numbers as vet 
We now seek to bring the second column of BA into agreement with the second column of 
In order not to change the form of the 


unit matrix by means of elementary transformations 
must avoid the following transformations 
b) interchanging the first and jth rows; 


Multiplying the first row by 


a 
adding @ times the 


st column, we 
onstant different from unity: ( 
st row to the jth row for a#0 

However, the remaining transformations are quite adequate to effect the desired reduction 


second column of BA is 
S07 


1067 


Confining attention to this column only, we see that the following sequence of transformations 


scalled for. We add a@ times the ith row to the jth row for the following succession of @, 7, and ) 


Performing these transformations on B gives C’, namely, 





two columns of CA are identical with those ofthe unit matrix The third colu 


11 


In dealing with this, we must 
we can reduce it to the desired form by adding @ times the 7th row to the jth row for the foll 


t now curtail transformations on the first two rows Neverthel 


Ing succession of a and j 


Performing these transformations on C gives D, namely 
1464 5481 6270 
350 1317 1509 
1373 S057 5745 


1447 


The first three columns of DA are identical with those of the unit matrix. The fourt! 


column is 
5 51061 


| 33644 
5 03483 
1 51308 
2058 
13 21462 
In dealing with this, we must now curtail transformations on the first three rows. Nevertheless 


we can start out to reduce it to the desired form by adding a times the ith row to the jth row for 


the following succession of a, 7, and j 





We now have the fourth column reduced to 


) 


annot reduce the 2 in the fifth row to unity except by dividing by 2, which would intro 
Iractions We would prete! to delay the introduction of fractions until the last possibli 


ent Accordingly, we multiply the first and third rows by 


then subtract the fifth row 
the first and third, then multiply the fifth row by 1. and, finally, interchange the fourth 
fifth rows 
We have now changed the first and third columns so that they are no longer the first and 
d columns of the unit matrix, but have a 2 where the unit matrix has a unit, Also oun 
rth column has the same property. This can readily be corrected by multiplving the appro 
e rows by one-half, but we will postpone this step until we can no longer avoid fractions 


Performing the indicated transformations on D gives the matrix 
17945 
HoH 
20337 
168720 


1554 


Chis is not exactly the matrix / as we defined it earlier, but is close enough so that we shall 


tk The first four columns of /-A are essentially the first four columns of the unit matrix 


having 2 in place of unity in the first, third, and fourth columns. The fifth column 


KA is 


347 
256 87577 


117 87645 


600 LO368 


We now reduce this bv adding a times the 7th row to the jth row for the following succession 


and j 





We now multiply the second, third, and fourth rows by 3, then add —1 times the sixth 
1 times the sixth row to the third and fourth. Finally, we intercha: 


row to the second and 


the fifth and sixth rows 
We could perform these transformations directly on EF to get a matrix that we will call 


However, it is computationally easier to proceed as follows. We perform the transformatio 


on the unit matrix, getting the matrix 


29 53877 5 80221 


108 36604 21 28600 
346 28689 68 02003 
339 7 66 75463 
90 19777 17 71726 
39 29215 


200 03456 


Then, if we multiply & on the left by the matrix above, we get the same matrix F' that wy 
would get by performing the indicated transformations on the rows of &. This matrix / 


follows 


$9935 13702 


1 83192 29419 


80075 96989 


2 93767 10487 


34504 23435 


1 26582 40724 


95298 44107 


3 49612 33901 


99845 26301 


3 66292 83279 


| 05671 28268 


3 87666 19778 


5 85396 37006 9 38741 43339 t 04497 80947 11 17196 52341 11 70499 53243 12 38798 74958 


— 5 74506 12736 9 21277 84365 3 96972 85793 10 96413 09887 11 48724 50163 


1 52478 89704 2 44515 11094 1 05360 03125 2 90997 52301 3 04881 42192 


3 38157 46308 5 42269 20056 2 33660 40527 6 45354 77403 6 76145 55311 7 15598 97507 


The matrix FA 
137 47502 92848 
504 34158 21303 
1611 63837 57000 
L581 65675 30712 
$19 78538 73419 


930 97185 49728 


Most of the transformations that we 


At this point, we can read off the determinant of A 
Compiling those that do 


used in forming F are such as to leave the determinant unchanged. 


we find that the determinant of F is 108. So the determinant of A is 


change the determinant 
5585 83112 98368 
One can readily write down the inverse of FA, namely, a matrix whose components ar 
fractions with the common denominator 930 97185 49728 and the following numerators 
165 48592 74864 0 0 68 73751 46424 


310 32395 16576 0 168 11386 0710! 
155 16197 58288 268 60639 59501 
155 16197 58288 263 60945 88452 
0 310 32395 16576 139 92846 24473 


0 0 1 





ally, we compute A~' from the equation A~'=(FA)"'F. This gives the same matrix for 
that was computed in section 1. 
The computations outlined in section 2 required 23 hours of computing time by a pro- 
sional computer. This time included the time needed to train the computer in the unfamiliar 
thod 


3. Remarks on Computational Details 


One advantage of the procedure outlined in section 2 is the ease in making numerical 
cks. Since all the operations are on rows, one can easily carry an extra check column 
ich is the sum of all the columns with which one is dealing. However, this is not needed, 
other checks are possible, as follows. The computation of B is easily checked by computing 
first column of BA and seeing if it agrees with the first column of the unit matrix. Then 
» can check C by computing the first two columns of CA and seeing if they agree with the 

st two columns of the unit matrix; and so on. 

Not only does this furnish a convenient check, but when a check is not forthcoming, one 
) often find the error by this method. For example, in computing FE, two mistakes were 

ide, and the resulting matrix was 


32262 1 47945 285 1 32919 13466 9968 


12391 HO699 % 4 54117 3115 38872 

3330 20337 317: 17997 1298 11972 

36052 1 68729 27507 51127 13242 

14752 74554 5727: 66587 2494 

83708 3 93352 2 975: 3 52225 29905 
hich differs from E in most elements of the fourth column (where the original error had snow- 
illed) and in the second element of the sixth column. When we multiplied this matrix on the 
cht by i, We got 


70652 
1 44778 
70650 
70650 
70650 


1 41300 


the first column, instead of twice the first column of the unit matrix. So, except for the 
ond row, our errors are all multiples of 70650, which is (450 157). As 157 is the fourth 
ment in the first column of A, it seems clear that there are errors in the fourth column of 
hat purports to be £, and that these errors are multiples of 450. With this information, the 
stake in computing in the fourth column was quickly discovered and eliminated. Now 
iltiplication on the right by A verified all but the second row, and a trivial amount of detec- 
work on the second row of the product sufficed to locate the error in the second row and 
its magnitude 
In setting up the computation for use on 1pm machinery, we notice that the majority 
the steps consist of adding @ times a row vector to another row vector, and that commonly 
sa small integer. It is not difficult to wire a multiplier so that if we insert a deck with a 
d containing a followed by cards with the components of the two vectors interleaved, the 
ltiplier will punch cards with the components of the resulting vector. At various stages 
the procedure, some columns are computed for the product of two matrices, of which the 


t has as rows the row vectors that we are manipulating, and of which the second is always A 
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Such a matrix multiplication can also be wired up for the multiplier. After the matrix m 


plication, one then makes a list of a’s and rows to be operated on by inspection from a colu 


This is most conveniently done by hand, which permits the exercise of judgment at this pe 


However, one can proceed perfectly mechanically, as we did in the illustrative example 


deed, it is not clear that one can really do much better by exercising judgment than we 
with our purely mechanical procedure Once the list of ae’s and rows to be operated oO 


compiled, the respective operations can be quickly performed on 1BnM machinery. W 


one must keep each row as a deck, and the row decks have to be interlea 


only a multiplier 
If a card programmed calculator is available, one may put 


and separated repeatedly 
and the operations 


entire row on a single card (unless the matrix is of reallv high order 


greatly speeded 


The matrix A that we used was constructed from a table of random numbers in an ef 


to furnish an example that might be considered typical 
With an increase in the order of the matrix to be inverted, the method presented he 
improvement over the standard methods. The method was 


becomes even more ot an 
Gunderson at Cornell Univers 


summer of 1948 while working with Dr. N. G 
in which we required the exact solution of 15 equations iu 
and a solution obtained 


devised in the 
on a proble m in number theory 


unknowns (one unknown was transposed to the right-hand side 
terms of it Fort inately, many of the coeflicients were zeros Even so, the usual meth 
of solution led us to hopelessly large numbers, whereas a solution was carried out by the met! 


) 


of this paper without encountering any integer of more than 12 digits 


Some ot the proced iral details of the present method were devised by DD: (csunders« 


The computations for the present paper were carried out by Lillian Forthal, Nancy Mar 


and Gerald Kimble, under the direction of Marvin Howard 


Los ANGELES st 14, 1950. 





